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Ibis invention is a device (10) for the percutaneous localization of a foreign object 
in a body has a base collar mount (16) and one or more movable collar mounts (12X14) all 
located on at least one perpendicular guide rail (18aX18b). The base collar mount (16) has 
a guide tube cannula (20) affixed to it. A series of tubing (28) and wire lengths (24) extend 
through the guide tube cannula (20) as the movable collar mounts (12X14) are advanced. As 
the tubing (28) and wire (24) are exposed from the guide tube cannula (20), they assume 
their native curvature thus forming a compound trajectory for the insertion of a foreign 
object in a body. 
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METHODS FOR STEREOTACTIC IMPLANTATION 



Government Support 
This work was supported by grant RR00166 
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Background Of The Invention 
The present invention provides methods for 
stereotactic implantation of foreign bodies in organs or 

15 tissues. More particularly, the present invention provides 

methods for stereotactic implantation of a foreign body along 
a compound straight and curvilinear trajectory through an 
organ or tissue. The compound trajectory may be identified by 
computer programs. The present invention allows short-term or 

20 long-term access to a disease locale by minimally invasive 
surgical techniques. 

Surgeons commonly treat localized disease processes 
by altering and, in some instances, by destroying biological 
tissue in a disease locale. In most cases, surgical treatment 

25 does not provide for chronic direct application of biological 
agents or other means of modulating or monitoring local 
tissues. Surgeries that involve such treatment have inherent 
limitations. In many clinical settings in which foreign 
matter is implanted for therapeutic purposes, such as hip 

30 replacements, major surgery is required. Overlying tissues 
and target tissues must be incised and opened to gain access 
to the implantation site. In some instances, surgeons may 
leave an implanted vascular catheter for localized 
administration of a drug. This type of treatment is generally 

35 limited, however, because the catheter delivers drug into a 

vess 1 that enhances delivery to th general circulation, the 
tumor may extend beyond the drainage of the v ssel, or 
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delivery is limited to the point of egress at the tip of the 
tubing. For example, localized chemotherapy for liver cancer 
can be provided by catheter izing the hepatic artery. The 
catheter must follow the route of hepatic artery. Infused 
5 drugs contact the entire drainage of the artery, eventually 
joining the general circulation in the inferior vena cava. 

In the treatment of brain lesions, such as 
malignancies, single vessels will often not perfuse the entire 
tumor bed. Further, the blood brain barrier acts to limit the 

10 delivery of drug across the capillary wall. As delivery of 

the drug is limited, the tumor receives low doses even during 
localized infusions of high doses. Attempts to directly apply 
medications to intracranial lesions are limited by the 
difficulties and dangers inherent in intracranial access. 

15 Localized radiotherapy of gliomas (brachytherapy) , for 

example, requires the passage of straight tubes through brain 
tissue to reach the tumor site. Often several such tubes are 
required to treat the entire tumor. Such treatments are 
limited because many delivery tubes may be required to reach 

20 the necessary points for insertion of the radioactive 

material. Also, linear percutaneous approaches to deep 
tissues are often impossible if the path traverses, or lies 
close to, critical structures that may be injured by catheter 
placement. 

25 Temporarily concealed retrieval hooks and loops of 

various designs have been used in angiography and urology. 
For example, O'Neill and Joseph describe retrieval of severed 
catheter tips from pulmonary arteries by means of a catheter 
from which a retractable loop could be extended (Affit HMrt Jti 

30 98, pp. 287-93, 1979). Similarly, ureteral stones have been 
snared and removed by catheters having extendable curved 
elements (Lassers, Am. Heart J. . 73, pp. 375, 1970). 

Looped and curved catheters of various designs have 
been used in angiography for placement of catheter tips into 

35 branched vessels. ( See , e.g. . Fellmeth et al., Radiology, 172, 
pp. 872-3, 1989; and Dorros and Levin, cathet. Cardiovasc. 
Diaan. . 13, pp. 138-140, 1987). These angiographic and 
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related techniques employed natural body channels for the 
placement of the catheters, i.e., blood vessels and the 
urinary tract. 

Concealed wire loops have been used in 
5 disconnection neurosurgeries, but for the inherent purpose of 
cutting tissue rather than penetrating with minimal damage as 
is desirable for treating localized lesions deep within normal 
parenchyma. The concealed curve is a loop that expands 
outward from a housing, thus shearing across tissue for the 

10 purpose of cutting fiber tracts (leukotomy) , rather than 

cleanly piercing along its trajectory. For example, such a 
technique has been employed for thalamotomy to relieve 
essential tremor (Blacker et al., J.Nerv. Ment. Pis.. 147, pp. 
49-55, 1968). Another type of probe has been introduced into 

15 the brain for a different purpose, without the intention of 
damaging diseased tissue. A flexible optical fiber tipped 
with a pressure transducer has been inserted into the brain 
for the measurement of intracranial pressure (Ostrup et al., 
J . Neuroaura. 67, pp. 206-207, 1987; Crutchfield et al., £l 

20 Neurosura. 72, pp. 482-487, 1990). Such devices have not been 
passed clear through the brain from the entry point to an exit 
point, and have not been used for drug application. No method 
is available for specifically selecting an extended region 
deep within the brain or body affected by a disease. 

25 Hany localised disease processes may be treated by 

systemically administered drugs. Oral or parenteral 
administration provides access to the general circulation and, 
thus, to the diseased site. Such systemic therapy is 
effective for certain diseases. Hany other diseases are not 

30 amenable to such treatment, however. Treatment may be limited 
by toxic side effects at other sites in the body. Treatment 
may be limited by the metabolism and elimination of the drug 
before it reaches the intended site. 

Molecular biology is leading to the emergence of 

35 new classes of therapeutic agents, many of which are 

bi logically too delicate and to xpensive to be useful as 
periph rally administered drugs. Means of delivering such 
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substances in clinically practicable ways are a major concern. 
This is esp cially true for the many currently untreatable 
diseases of the brain. Therapeutics often cannot cross the 
blood brain barrier in sufficient quantities to produce 
5 effective tissue concentrations. Likewise from a research 

perspective, the understanding of brain function is impeded by 
the fact that substances can be administered locally only to a 
single "point site", and only acutely or subchronically. 
Long-term studies of localized drug effects are hampered by 

10 the inflammatory processes associated with normal implanted 
cannulae, which are stiff and usually metallic. The 
experimental implantation of encapsulated tumor cells that 
produce neuroactive substances (Erne rich et al., Weurosci . 
Biobehav. Rev. . 16 , pp. 437-447, 1992) is limited by inherent 

15 dangers of implanting tumor cells, by the lack of means to 
regulate or to alter the flow of substances after 
implantation, and by the limitation to point sites of 
implantation. The effects of drugs in anatomically distant 
brain locales are not easily tested by systemic delivery or by 

20 existing means of administration to a point site. 

What is needed in the art are means for accessing 
extended local regions of organs, which may be vascularized by 
several different and often inaccessible vessels. Preferably, 
the means would allow for accurate percutaneous access to 

25 disease sites that would not compromise sensitive structures, 
such as points in the brain stem, as well as allow for chronic 
localized therapy of a lesion. Quite surprisingly, the 
present invention fulfills these and other needs. 

30 ffMTfflfimrY Qf The Invention 

The present invention provides devices having a 
proximal end and a distal end for percutaneous localization of 
a foreign object in a body (Figs. 1A-1C and 2A-2E) . The 
devices generally comprise a plate-like base collar mount 

35 having a proximal surface, a distal surface rigidly affixed to 
at least on perpendicular guide rail, and a chann 1 between 
the proximal and distal surfaces through which channel a 
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guide-tube cannula extends from the proximal surfac of the 
base collar; a plate-like first movable collar mount slidably 
located on each guide rail rigidly affixed to a curved first 
extendable cannula having a proximal end and a distal end, 
wherein the distal end of the first extendable cannula is 
affixed to the first movable collar in a channel through the 
first movable collar and the proximal portion of the first 
extendable cannula is in the guide-tube cannula; and a plate- 
like second movable collar mount slidably located on each 
guide rail rigidly affixed to a second extendable cannula or 
wire having a proximal end and a distal end, wherein the 
distal end of the second extendable cannula or wire is affixed 
to the second movable collar in a channel through the second 
movable collar and the proximal portion of the second 
extendable cannula or wire is in the first extendable cannula; 
wherein proximal sliding of the first movable collar will 
extend the first extendable cannula through the proximal end 
of the guide-tube cannula and proximal sliding of the second 
movable collar will extend the second extendable cannula or 
wire through the proximal end of the first extendable cannula* 
The devices may have additional movable collar mounts and 
extendable cannulae. Generally, the path of percutaneous 
access will be a compound trajectory comprising straight and 
curvilinear segments. 

Methods for percutaneous placement of foreign 
substances at desired locations in the body are also provided. 
The methods generally comprise positioning a device as 
described above at a predetermined location and orientation 
over the body; advancing the guide-tube cannula into the body, 
covering an original predetermined distance; advancing the 
first extendable cannula into the body, covering a first 
additional predetermined distance; optionally advancing a 
second or more additional coaxial cannulae, covering a second 
or more additional predetermined distances; and introducing 
the foreign substance through the final extendable cannula to 
the desired location in the body. Th d sired position in the 
body is gen rally id ntif led by a diagnostic imaging techniqu 
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(e.g., magnetic resonance imaging, computerized tomography, 
ultrasound, or th like) . The foreign substanc is often a 
pharmaceutical compound or catheter. 

Methods for treating localized diseased regions in 
5 a body are also provided. The methods involve 

stereotactically implanting a foreign object along a compound 
trajectory through a desired region in a host. The methods 
generally comprise identifying an arcuate or helical 
curvilinear path through the host, which curvilinear path may 

10 traverse the region to be treated; inserting a guide-tube 

cannula through a surface of the host to meet the curvilinear 
path; advancing a curved cannula or wire, preformed to match 
the curvilinear path, through the guide-tube cannula into the 
host along the curvilinear path to a predetermined point; 

15 optionally advancing one or more additional coaxial cannulae, 
the last of which may be a solid wire, through the curved 
cannula to the surface of the host; attaching a foreign object 
to the final cannula or wire; withdrawing each wire or cannula 
through the host in reverse order, whereby the foreign object 

20 traverses the compound trajectory including the insertion path 
of each wire or cannula in the host. The foreign object may 
include a cryogenic probe, an electrode or electric device 
such as chips or wires, or a semipermeable catheter. The 
catheter may be used for analyzing substances passing into the 

25 catheter from the brain, for exposing the brain region to 

therapeutic compounds elaborated by cells introduced into the 
catheter, or for adjustably exposing the brain region to 
therapeutic compounds placed directly in the catheter itself. 
The treated disease may be a tumor, Parkinson's disease, 

30 schizophrenia, depression, Alzheimer's disease, epilepsy, 

spinal cord trauma, or the like. The therapeutic compound may 
be a dopamine agonist, a neurotrophic factor, or any of a wide 
variety of other agents. 

Kits for placing catheters along compound 

35 trajectories in a host are also provided. The kits generally 
comprise a d vie of th pr sent inv nti n and a cath ter. 
Th catheter is typically semiperm able. Often the catheter 
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is polyacrylonitrile polyvinylchloride copolymer* Pumps may 
also be supplied in the kits. 

Methods for determining a compound trajectory 
through a localized or extended region of a body for selective 
5 placement of a therapeutic modality are also provided. The 
methods generally comprise identifying the body region in a 
diagnostic image of tissue surrounding the region; assigning a 
coordinate system to the body; identifying three or more non- 
coil inear points in the region; identifying the circle that 

10 intersects the three points, or a helix that passes near the 
four or more points; determining the radius, center and plane 
of the circle, or the formula of the helix; identifying points 
of intersection of the circle or helix with predetermined 
parallel sections through the body; identifying two points of 

15 tangency at which entry and exit segments of the compound 

trajectory contact an arc of the circle or helix; determining 
the direction cosines of the entry and exit segments; 
calculating the arcuate distance between each pair of circle 
or helix intersection points on consecutive sections, from one 

20 point of tangency to the other; summing the arcuate distances; 
identifying points of intersection of the entry and exit 
segments with the predetermined parallel sections through the 
body; calculating the linear distances between each pair of 
entry or exit segment intersection points on consecutive 

25 sections; summing the linear distances; determining the 

compound trajectory of the path as the points of intersection 
of the arcuate and linear segments with the predetermined 
parallel sections; redefining the trajectory in the coordinate 
system of the diagnostic image; and displaying the trajectory 

30 in the diagnostic image. The diagnostic image is generally a 
magnetic resonance image, computerized tomographic image, 
ultrasound image, or the like. The body region is often in 
the brain, and may be, for example, a brain tumor or a portion 
of the basal ganglia. 



35 
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p r ief Description of the Drawings 
Figs* 1A - 1C illustrate a d vie f th present 
invention having two cannulae and one wire, in the fully 
extended position, partially extended position, and non~ 
5 extended position. 

Figs. 2A - 2E illustrate a device of the present 
invention having four cannulae and one wire, in the fully 
extended position, partially extended positions, and non- 
extended position. 
10 Fig. 3 illustrates the position of a catheter 

placed by the methods of the present invention for treatment 
of a deep site in the brain. 

Fig. 4 illustrates a perspective view of a device 
of the present invention having two cannulae and one wire in 
15 the fully extended position, mounted on a Kopf stereotactic 

instrument with PED and ROT adjustments at the standard square 
and upright starting positions. 

Fig. 5 illustrates a partial cut-away view of the 
device of Fig. 1, emphasizing the initial curvature of the 
20 curved cannula within the straight guide-tube cannula, and the 
initial straight portion of the innermost wire within the 
curved cannula (Description of Specific Embodiment, Part V, 
Step 1) . 

Fig. 6A-B illustrate geometrically the PED (6A) and 
25 ROT (6B) adjustments for a device of the present invention, as 
mounted on a Kopf stereotactic instrument (Description of 
Specific Embodiment, Part V, Step 3). 

Fig. 7 illustrates geometrically the AX adjustment 
along the insertion axis for a device of the present 
30 invention, as mounted on a Kopf stereotactic instrument 
(Description of Specific Embodiment, Part VI, Step 1). 

Fig. 8 illustrates geometrically the CR adjustment 
across the insertion axis for a device of the present 
invention, as mounted on a Kopf stereotactic instrument 
35 (Description of Specific Embodiment, Part VI, Step 2). 

Fig. 9 illustrates geometrically the calculation of 
the length of the arc al ng the curvilinear portion f the 
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trajectory from one inters cting parallel plane to the next, 
for an arcuate curvilinear trajectory, without helical torsion 
(Description of Specific Embodiment, Part VII, Step 1). 

Pig. 10 illustrates geometrically the calculation 
5 of the most posterior and anterior points on the circle in 
which the curvilinear portion of the trajectory resides, in 
the case of an arcuate curvilinear trajectory, without helical 
torsion (Description of Specific Embodiment, Part VII, Step 
4). 

10 Fig. 11A-B depict quantitative records of turns 

made contralateral to the side of a lesion and an implant, in 
response to drug treatment via the implant. The number of 
turns (Y-axis) completed each minute throughout the 24 hours 
of the day (X-axis) is indicated by a line of appropriate 

15 height. 

Descriptio n Of The ffp?fflfiff Embodiments 
The present invention provides devices and methods 
for determining curvilinear paths through three dimensional 

20 structures. The methods allow clinicians to select a targeted 
group of points in a compact or extended region identified in 
an imaging study, such as MRI, CT, or ultrasound, and to link 
the points by a helical curvilinear path. The targeted region 
may be remote from the surface of the body. The region may 

25 have a tortuous shape, or may need treatment near its margins, 
requiring access from a helical arc or a coil of several 
periods. In other cases the region will have a simpler shape 
and will require access centrally, so that often the helix 
will be circular. After determining the helical path, devices 

30 of the present invention may be constructed to provide a means 
for implanting a catheter through the targeted region and 
treating the targeted region including each of the selected 
points. The parameters of the helix provide the parameters 
employed in shaping the curved cannulae as described above. 

35 One specific embodiment of the present invention 

will now be described by ref rring t Pig. 1, where an 
implantation tool 10 is shown. Th implantation tool 10, in 
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one embodiment:, includes an upper collar mount 12, a middle 
collar mount 14, and a low r collar mount 16. Th se collar 
mounts are each formed, e.g., out of plastic and are plate- 
like in shape. The mounts serve as handles and as guides for 
5 mounting and manipulating other components of the implantation 
tool 10. The lover collar mount 16 is fixed to a pair of 
stainless steel guide rails I8a-b, while the upper and middle 
collar mounts 12, 14 are slidably connected to the guide 
rails. 

10 The tool further consists of a guide-tube cannula 

20 which is affixed to the lower collar mount 16 and which 
slidably receives a length of curved nitinol tubing 28. 
Although other materials may be used for the tubing and wiring 
of the tool, nitinol, a nickle-titanium alloy, is preferred as 

15 it exhibits gxeat flexibility and a strong shape memory. The 
length of curved nitinol tubing 28 is, in turn, affixed to the 
middle collar mount 14 and slidably receives a length of 
nitinol wire 24. The length of nitinol wire 24 is attached to 
the upper collar mount 12. A gauge 26 for the curved nitinol 

20 tubing 28 is provided on the upper collar mount 12, slidably 
fitting around the guide rail 18a. Another gauge 22 is 
provided for the nitinol wire 24. The gauge 22 is affixed to 
the middle collar mount 14 and adjustably affixed to the upper 
collar mount 12. Finally, a stop 30, attached to the middle 

25 collar mount 14, is used for the curved nitinol tubing 28. 

The basic operation of an implantation tool 10 is 
shown in Figs. 1A-C. In Fig. 1A the tool 10 is positioned for 
the mflyimw extension of the nitinol wire 24 and tubing 28. 
As the upper collar mount 12 moves away from the middle collar 

30 mount 14, as in Fig. IB, the nitinol wire 24 retracts into the 
curved nitinol tubing 28. When the upper collar mount 12 
reaches the end of the gauge 22, it is tightened in position 
on the gauge 22. When the upper collar mount 12 and the 
middle collar mount 14 move together away from the lower 

35 collar mount 16, the curved nitinol tubing 28 and the nitinol 
wire 24 retract into the guide- tub cannula 20. The gauges 
nsure that the tubing and wir do not retract to far. The 
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gauge 22 also ensures that the upp r and middle collar mounts 
12, 14 move in unison when the curved nitinol tubing 28 and 
the nitinol wire 24 are retracted together into the guide- tube 
cannula 20. The properties of nitinol allow the curvature of 
5 the curved tubing and the wire to be precisely controlled, 

even after residing within the straight guide-tube cannula 20. 
As the tool 10 is extended into the body, the tubing and wire 
emerge from the guide-tube cannula 20 and express their native 
curvature and the device thus defines a distinctive pre- 

10 planned trajectory. The curvature may be programmed into the 
nitinol using means well-known in the art. 

At maximal extension, the nitinol wire protrudes 
from the body. Semipermeable tubing can be fastened to its 
tip so that when the tool 10 is retracted the semipermeable 

15 tubing is drawn through the body along the pre-planned 
trajectory. 

Those skilled in the art will appreciate that the 
curvatures and relative lengths of the wires and tubing may be 
modified as needed for specific applications. Skilled 

20 practitioners will also recognize that the implantation tool 
10 may be used in conjunction with a variety of stereotactic 
instruments, such as a Kopf stereotactic device as shown in 
Fig. 4. The Kopf stereotactic device is available from David 
Kopf Instruments, Tujunga California, and is used as an 

25 example platform in the discussion herein. 

Another specific embodiment of an implantation tool 
40 is shown and described by referring to Figs. 2A-E. In this 
embodiment, five elements are utilized to provide additional 
trajectory options for reaching remote and otherwise 

30 inoperable sites in the brain or other parts of the body. The 
tool includes five collar mounts 42-50. The top four collar 
mounts 42-48 are slidably attached to two stainless steel 
guide rails 41a-b. The bottom collar mount 50 is fixed to the 
ends of the guide rails 41a-b. A guide-tube cannula 52 is 

35 attached to the bottom collar mount 50 and slidably receives a 
first 1 ngth of curved nitin 1 tubing 56. This first length 
f curved tubing 56 is attached to th n xt highest collar 
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mount 48 and receives a length of straight nitinol tubing 58. 
The straight tubing is attached to the next highest collar 
mount 46 and receives a second length of curved nitinol tubing 
60 which in turn is attached to the next highest collar mount 
5 44. A length of nitinol wire 54 slides within the second 

length of curved nitinol tubing 60. The length of wire 54 is 
attached to the. upper collar mount 42. Thus, as is depicted 
in Pig. 2A, as the collar mounts are moved together, the wire 
54, the second curved nitinol tubing 60, the straight nitinol 
10 tubing 58, and the first curved nitinol tubing 56 are all 

extended and exposed past the end of the guide-tube cannula 
52. When the collar mounts are moved apart, as shown in Fig. 
2E, the tubing and wire are retracted within the guide-tube 
cannula . 

15 The implantation tool 40 also includes a variety of 

gauges attached to each of the collar mounts. For example, 
the length of straight nitinol tubing 58, which is attached to 
collar mount 46, has a corresponding gauge 66, which is 
adjustably attached to collar mount 46. Further, the gauge is 
20 also attached to the next-lowest collar mount 48. Observation 
of the gauge 66 as the collar mount 46 is raised enables the 
length of straight nitinol tubing 58 to be retracted the 
proper distance. Its tip is flush with the tip of the first 
curved nitinol tubing 56 when the top of its collar mount 46 
25 is flush with the top of the gauge 66, as shown in Fig. 2D. 
Tightening the attachment of collar mount 46 on gauge 66 
ensures that collar mounts 46 and 48 will subsequently move in 
unison during the retraction of the first curved nitinol 
tubing 56. Each of the other lengths of nitinol has a 
30 corresponding gauge. The nitinol wire 54 is gauged by gauge 

62. The first and second lengths of curved nitinol tubing are 
gauged by gauges 68 and 64, respectively. Those skilled in 
the art will realize that the lengths of the gauges will be 
chosen to correspond to the necessary lengths of nitinol 
35 tubing and wiring required for a specific procedure. The 

calculation and selection f such lengths and curvatures will 
be discussed infra in conjunction with a description of the 
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methods for planning a trajectory of the thre part 
implantati n tool. 

Referring now to Fig. 3, a view of a human cranium 
having an implant inserted by a device according to the 
5 present invention is shown. The figure relates the compound 
trajectory which must be followed in order to make such an 
implant. An implant with a front, non-permeable portion is 
shown as item 90. The implant also includes a back non- 
permeable portion 92 and a semipermeable portion 94. The 

10 interrupted line format of portions 90 , 92, 94 in the drawing 
indicate that they pass deeply through the substance of the 
brain, rather than lying at its surface. Return and input 
lines 96, 98 are also shown which connect the example implant 
with a pump. These lines are solid in the drawing, to 

15 indicate that they lie outside the brain, just under the skin. 
Those skilled in the art will appreciate that other foreign 
objects may also be inserted using the device and method of 
the present invention, including a flexible catheter ending 
deep in the brain or bbdy without a return line, various 

20 pharmaceutical compounds, cryogenic probes, or the like. 

For example, in one embodiment the foreign object 
is a compound catheter. A central portion of the compound 
catheter is attached to connecting portions, which are 
inserted smoothly inside the central portion to a preselected 

25 point. A small connector bead (e.g., 25-thin stainless steel 
tubing cut into smooth lengths) , optionally coated with 
Teflon m or the like, is lodged in the inserted tip of each 
connecting portion. The bead has a diameter large enough to 
press the outer wall of the connecting portion against the 

30 inner wall of the central portion, thereby producing a tight, 
friction-fit seal. The friction-fit seal is produced using, 
for example, a fine thread or suture, passed through the 
central and connecting portions and the beads, tied in a small 
slip-knot between the two beads, and pulled in succession from 

35 each end against the corresponding bead to bring it into the 
tip of th conn cting portion, then pulled at both ends to 
remove th slip-knot and subsequently the sutur • A small 



WO 96/03087 



PCTYUS95/09391 



14 

terminal bead is also friction-fit into one end of the 
compound catheter. The terminal bead preferably has each end 
internally shaped to a serrated conical form, allowing a firm 
friction fit to a complementarily shaped conical tip of a 
5 wire. A fine wire (e.g., 4 -mil stainless steel) is then 

passed through the entire length of the compound catheter and 
lodged at its tip in the terminal bead. This assemblage 
provides a stronger implant. 

In an another aspect the method and apparatus for 

10 stereotactically implanting a foreign object along a compound 
trajectory in a brain or body further comprises a coordinated 
means for pushing the implant into the brain or body from 
above as it is pulled from below. This aspect involves the 
coordinated withdrawal of each wire or cannula in reverse 

15 order. Two concentric implant-bearing tubes of differing 

diameter and capable of telescoping with the other are used to 
encase the implant. The tube having the larger diameter is 
fastened to a block containing a pulley mounted on a 
stereotactic carrier, so that the tube can be lowered over the 

20 smaller tube. Thus, the encased implant can be advanced, 
under the guidance of the larger tube, into the insertion 
trajectory in the brain or body by the planned distance. The 
pulley-wheel is fastened to the top of the stereotactic 
carrier on a rotational ly adjustable mount, such that its axis 

25 is perpendicular to the line of the two concentric telescoping 
tubes and its plane of revolution can be adjusted to include 
the top of a second pulley-wheel on a rotationally adjustable 
mount fastened to the top of the device bearing the guide-tube 
cannula, curved first extendable cannula, and second 

30 extendable cannula or wire for insertion. A non-stretchable 
pulley cable or wire is fastened at one end to the block that 
drives the larger implant-bearing tube, and runs up from the 
block to the pulley-wheel fastened to the top of the 
stereotactic carrier. The flexible wire then extends across 

35 to the pulley-wheel at the top of the stereotactic carrier 

bearing th guide-tube cannula, curved cannula, and wire, and 
finally xtends down parallel t the axis of the guide-tube 
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cannula to a fixation point in the second or outermost collar- 
mount which bears the innermost wire. A solid brace of 
adjustable length extends between the mounts of the 
rotationally adjustable pulley-wheels to maintain their 
5 positions when the non-stretchable pulley cable or wire is 
pulled. The motion of the larger implant-bearing tube 
simultaneously and at the same rate lowers the implant into 
the brain or body from above and raises the collar-mount which 
bears the innermost wire, thereby pulling the implant into the 

10 brain with the same motion. 

Methods for planning a trajectory for a device of 
the present invention are described below. The implantation 
tool is temporarily advanced in stages along the planned 
compound trajectory, into the brain or another organ or region 

15 of the body, until its leading tip emerges again; then a 

length of implantable tubing is connected to the tip and drawn 
into position as the telescoping implantation tool is 
retracted. Various software packages were utilized in 
conjunction with a portion of the steps to be described. 

20 Those skilled in the art will recognize that methods for the 
calculation and planning of such a trajectory may differ in 
detail from those utilized. 

The methods will be described in nine parts. The 
first part (Machine-dependent Input) will begin with imaging 

25 scans, proceed to three original target points, and add other 
measurements and user-input adjustments to provide to the Core 
Section. The next two parts (Core Section) will begin with 
the three original target points, incorporate two additional 
user- input tangent points, and proceed to the final series of 

30 calculated points all along the trajectory. The fourth, 

fifth, sixth, and seventh parts. (Machine-dependent Output) 
will begin with the calculated trajectory and proceed to its 
display, its improvement and recalculation if necessary, the 
angular and translational settings it prescribes for the 

35 stereotactic equipment, and the specifications it provides for 
th implantabl tubing and the implantation tool. Th ighth 
part (Alternative Input Methods) will return to the beginning 
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and describe a few of the many possible alternative input 
techniques that have b en developed in the prototype. The 
ninth part (Fitting a Three-Dimensional Helical Curve to a Set 
of Points) will present an alternative to the second part, for 
5 the case in which a three-dimensional helix is desired and 
more than three target points have been entered* 

I. Machine-dependent input Algorithms: From an imaging 
scanner (MRI, CAT, ultrasound, etc.), to an accessible file 
10 format representing its images, to a display with a graphical 
user interface (GUI) , to inputs in the form of selected 
pixels, to the selection of three original target points. 
Additional hardware parameters and user- input adjustments. 

15 1. A first step is to establish physically an 

orthogonal coordinate system. The manner by which this system 
is established is immaterial, as long as it can be effectively 
reproduced in both the radiological and surgical environments. 
It is necessary to configure equipment that can be used to 

20 align the imaging and surgical coordinate systems. To 

accomplish an accurate placement of implanted tubing, it is 
important for the body to be in the same fixed position in the 
imaging machine and in the surgical context. There are a 
number of ways to align the two systems, depending on the 

25 targeted part of the body and on the appropriate and available 
equipment. One useful way is illustrated below. 

The example described herein reflects the 
development of the technique of the present invention in the 
context of nonhuman primate neurosurgery, making use of a 

30 standard stereotactic instrument fixed at the external 

auditory meati and the inferior orbital ridge. Those skilled 
in the art will appreciate that the method and apparatus of 
the present invention is preferably used in conjunction with 
tools allowing implantation in human organs and tissue. The 

35 nonhuman primate stereotactic instrument is described as a 
c nvenient example. 
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The origin at which this system is zeroed (called 
earbar zero) is a point midway between the points at which the 
stereotactic earbars contact the bony external auditory 
meatus. The line between these two contact points is one of 
5 the axes (the ML^axis) of the coordinate system. A second 
axis (the AP axis) passes through earbar zero and a point on 
the sagittal pla ne at the level of the inferior orbital ridge 
(its superior surface). Together, the ML axis and the AP axis 
define a standard plane, which is known to classical 
10 craniometry as the Frankfurt plane, or to tomographic ^ 

radiology as an axial plane. The third axis (the DV axis) is 
orthogonal to both the ML and AP axes, passing through earbar 
zero and through the cranium above earbar zero* The plane 
orthogonal to the ML axis is the midline sagittal plane. The 
15 plane orthogonal to the AP axis is known to tomographic 
radiology as the horizontal plane , or to classical 
neuroanatomy as the coronal plane. The utility of the method 
and the algorithms is not limited, however, to the brain 
stereotactic context. More generally, in the summary 
20 statements and in the formula derivations presented below, the 
following terminology conventions are respected: 
earbar zero » origin 
ML axis « X axis 
DV axis - Y axis 
25 AP axis - Z axis 

coronal plane - DV/ML plane « XY plane 
sagittal plane - DV/AP plane * Y2 plane 
axial plane « ML/AP plane » XZ plane 

30 The following sign conventions are also respected 

in the formula derivations presented below. The X values 
along the ML axis are positive to the left of earbar zero, 
negative to the right. The Y values along the DV axis are 
positive above earbar zero, negative below. The Z values 

35 along the AP axis are positive in front of earbar zero, 
negative behind. The c nvention that requires this sign 
designation was set out in rder to c nf rm to the standard 
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examination of coronal sections in neuroanatomy. It permits 
the surgeon to examine coronal sections as if viewing the 
brain from in front of the face, and in that situation to have 
the view on the right (patient 1 s left) positive and the view 
5 on the left (patient 1 s right) negatiye. 

The stereotactic device (Kopf stereotactic 
instrument, #1404, David Kopf Instruments, Tujunga CA) (Fig. 
4) fixes the cranium in place by four contacts, two at the 
bony external auditory meati (earbars) and two fixed against 

10 the inferior orbital ridges of the eyes (eyebars) . The 

stereotactic device for nonhuman primate magnetic resonance 
imaging (MRI) developed at the Regional Primate Research 
Center (RPRC, University of Washington, Seattle, WA) fixes the 
cranium in the same manner. In both cases, the plane in which 

15 these four bony landmarks reside (the Frankfurt plane) is 
typically held parallel to the plane of the arms of the 
stereotactic instrument, but it can be adjusted if necessary 
by rotating the cranium on the earbars and adjusting the 
eyebars (Dubach et al., J- Neurosci, Meth. . 13, pp. 163-169, 

20 1985) to a higher or lower position. As a rule, the eyebars 
are set to the same height as the earbars, but other 
variations are acceptable as long as the settings can be 
reproduced on both devices; the bony landmarks must be in the 
same position for scanning and for surgery. The stereotactic 

25 device for the imaging scanner is mounted in the imaging field 
in a fixed position such that, if eyebars and earbars are set 
to the same height, the Frankfurt plane containing these 
landmarks is parallel to the axial plane of the imaging field. 
The mount is also designed to place the cranium's midline 

30 sagittal plane in the center of the imaging field. Therefore, 
if the cranium were properly and perfectly positioned, and if 
the imaging scanner always centered its images at the same 
point, no registration points would need to be entered. 

35 2. A second step is to configure equipment that 

can be used to interfac betwe n the imaging machine and the 
analytic algorithms that plan th trajectory. Again, th re 
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are a number of ways for accomplishing this. One us ful way 
is illustrated below. 

When the imaging machine generates a video picture, 
that picture is the expression of an underlying data file. 
5 Such a file, converted to a TIF file in the rectangular format 
as described below, is a series of numbers (of magnitude 0-254 
inclusive) representing gray-scale values of pixels. The 
first pixel in the data portion of the file is the pixel in 
the upper left-hand corner, the second is adjacent to its 

10 right, the third is to the right of that, and so on to the end 
of a line of 256 pixels. The 257th pixel is on the left just 
below the first pixel, starting another row of 256 pixels 
proceeding to the right. There are 256 such rows, from the 
top to the bottom. The origins of the imaging stereotactic 

15 coordinate system and the surgical stereotactic coordinate 

system, on the other hand, are in the center. Because of the 
structure of a TIF file, it is most straightforward to define 
for the scan a TIF coordinate system with its origin at the 
upper left-hand corner of the picture, with X values (ML) 

20 increasing toward the right of the picture, and Y values (DV) 
increasing toward the bottom of the picture. 

In one specific embodiment, software may be 
utilized in order to assist in the selection of input points. 
For example, the following may be used: 

25 (1) A conversion program available as part of the 

imaging facility converts each file, containing data from a 
coronal tomographic slice, from its native imaging scanner 
format to a rectangular TIF file format. 

(2) An FTP network protocol, run from an 80486 

30 microcomputer r may transfer the file onto the 80486 computer. 

(3) An 80486 conversion program (such as those 
available at the RPRC, implemented in Microsoft Visual BASIC) 
can be used to put the data into a 256x256 format, and to 
query the user to set the gray-scale in such a way as to make 

35 full use of 255 shades of gray (omitting use of the brightest 
shade) . 
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(4) The file nay then be displayed by a commercial 
image-display package (Adobe Photoshop) , for visual 
observation* Any pixel in the grid can be highlighted to the 
brightest gray-scale shade. 
5 (5) An 80486 program (which may, for example, be 

implemented in Microsoft Visual BASIC) can then be used to 
survey the modified TIF file to find highlighted points, 
convert each point from its sequential position in the TIF 
file to its TIF coordinates as defined above, and pass the 
10 list of points into a text file- The text file may then be 
utilized as input information for analysis. 

3. A third step, given these configurations, is to 
register the imaging stereotactic coordinate system with the 

15 surgical stereotactic coordinate system assumed by the 

subsequent algorithms. This registration is accomplished by 
informing the algorithms of the location of landmarks in the 
imaging system, so that input from and output to images can be 
properly expressed within the algorithms. 

20 In practice, the imaging machine does not center 

its images at a standard single point, so at least one point 
is necessary for registration. In practice, furthermore, the 
positioning of the subject is not perfect, so that more points 
are required for registration. A total of eight registration 

25 points are entered. An imaging stereotactic device may 

contain channels of oil in the earbars, which can be perceived 
in the pictures it creates. On the basis of a display on the 
computer monitor, a single point (LeftEarML, LeftEarDV, 
LeftEarAP) is entered for the left earbar, and another 

30 (RightEarML, RightEarDV, RightEarAP) for the right earbar, 

marking the coronal slice and the ML/DV location of the pixel 
nearest to the corner of the earbar, which is more specifiable 
and more faithfully rendered by the imaging scanner than other 
points along the channels. It is preferred, however, to 

35 determine a set of landmarks from a position close to the 
targeted p sit ion; th earbar points are well utside the 
cranium, and th location of th midline sagittal plane in th 



WO 96/03087 



PCTAJS95/0939J 



21 

image can vary significantly from one recording session to 
another. Therefore, tw p ints ar ntered al ng the midline 
of the coronal section through each original target point, one 
point at a relatively dorsal position along the midline (e.g., 
5 UpMidPlML, UpMidPlDV, UpMidPlAP) and the other at a relatively 
ventral position along the midline (e.g., LoMidPlML, 
LoMidPlDV, LoMidPlAP) • 

4. A fourth step is to enter the three original 
10 target points. 

There are many different ways to proceed in this 
step. (See Part VIII for some alternatives.) One way is to 
acquire the three original target points, (PT1ML, PT1DV, 
PT1AP) , (PT2ML, PT2DV, PT2AP) , and (PT3ML, PT3DV, PT3AP) , in 
15 the manner described above for. the earbar points. Along with 
each target point, two points on the midline in each 
corresponding section are also acquired, as described above. 

The points acquired in Steps 3 and 4 can be 
summarized: 

20 (LeftEarML, LeftEarDV, LeftEarAP) 

(RightEarML, RightEarDV, RightEarAP) 



(PT1ML, PT1DV, PT1AP) 
(UpMidPlML, UpMidPlDV, UpMidPlAP) 
25 (LoMidPlML, LoMidPlDV, LoMidPlAP) 

(PT2ML, PT2DV, PT2AP) 
(UpMidP2ML, UpMidP2DV, UpMidP2AP) 
(LoMidP2ML, LoMidP2DV, LoMidP2AP) 

30 

(PT3ML, PT3DV, PT3AP) 
(UpMidP3ML, UpMidP3DV, UpMidP3AP) 
(LoMidP3ML, LoMidP3DV, LoMidP3AP) 



35 



5. A fifth step is to convert the input points 
from pix 1 wits to millimet r units and to place them in th 
ster otactic coordinate system used by the Core Algorithms. 
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Again, there are a numb r f ways to do this. On 
useful way is illustrated b low. 

The tilt and translation of the imaging coordinate 
system is detected and compensated on the basis of the line 
5 between the eyebars and the midline of each slice containing 
an original target point. First, for the earbars, in the TIF 
coordinate system, the line (in the coronal plane) connecting 
the two earbar points has the general formula: 

10 (1) 

y - m * X + b 

or: 
(2) 

15 y - BarbarSlope * X + Earbarlntercept 

The slope of the line is calculated as (aY / aX) , 
the change in Y from one entered earbar point to the other, 
divided by the change in X: 

(3) 

BarbarSlope « (LeftBarDV - RightBarDV) / (LeftEarML - RightBarML) 



20 



The intercept (Earbar Intercept) with the Y axis 
25 (the left-hand margin of the TIF file, off to the right of the 
brain) can be derived by substituting (3) in (2): 

(4) 

Y - ((LeftBarDV - RightBarDV) / (LeftEarML - RightBarML)) 
30 * X + Earbarlntercept 

rearranging to obtain: 
(5) 

35 Earbarlntercept - Y - ((LeftBarDV - RightBarDV) / 

(LeftEarML - RightBarML)) * X 
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and then substituting X and Y valu s from one of the earbar 
points (LeftEarML, LeftEarDV, LeftEarAP) in (5): 

(6) 

Earbarlntercept - Lef tBarDV - ( (LeftEarDV - RightEarDV) / 
(LeftBarML - RightBarML)) * LeftEarML 

The same derivation provides the coefficients in 
the formula for the line connecting the upper and lover points 
entered from the midline on the slice containing the second 
original target point (or, in the same manner, the midlines of 
the slices containing the first and third original target 
points) . 

(7) 

MidP2Slope - (UpMidP2DV - LoM±dP2DV) / (UpMidP2ML - LoMidP2ML) 

(8) 

MidP2 Intercept - LoMidP2DV - <<UpMidP2DV - LoMidP2DV) / 
(UpMidP2ML - LoMidP2ML) ) * LoMidP2ML 

The AP position of the earbars in the MR! picture 
is calculated as the average of the two earbar AP values: 

(9) 

BBZeroAP * RightBarAP + (LeftBarAP - RightBarAP) / 2 

The field of view (FOV) in centimeters and the 
slice thickness (SliceThickness) of the tomographic slices are 
provided by the imaging scanner* Given that the standard 
display in the TIF file is a 256x256 grid, and that it 
represents the entire field of view (FOV x FOV centimeters) , 
the edge of each square pixel in the grid can be calculated in 
millimeters: 

(10) 

PixelSiae « (FOV * 10) / 256 
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The position of th s cond targ t point can now be 
translated into terms of millimeters from earbar zero. Its AP 
position is measured as the distance of its section from 
EBZeroAP, and is adjusted by the SliceThickness factor: 

5 

(11) 

AP B - SliceThickness * (Pt2AP - EBZeroAP) 

Its ML position within its section is defined as 
10 its distance to the midline, the line between the points 

(UpMidP2ML, UpMidP2DV, UpMidP2AP) and (LoMidP2ML, LoMidP2DV, 
LoMidP2AP) . This distance can be written in general , for the 
point (x 0 , y 0 ) and the line (AX + BY + C = 0) as: 

15 (12) 

d - | Ax Q + By Q + C | / SQRT(A 2 ♦ B 2 ) 
(See, Taylor, Calculus With Analytic Geometry- Prentice Hall, 
Englewood Cliffs, NJ, p. 272, 1959) 



20 The formula for a line used above in the derivation 

of the slope and intercept of the midline, rearranged from 
(1), was: 

(13) 

25 m*X-Y + b»0 

In this formula as written, the coefficients (from 
(7) and (8)) are: 



30 (14) 

A - MidP2Slope - <UpMidP2DV - LoMidP2DV) / (UpMidP2ML - LoMidP2ML) 

(15) 

B - -1 

35 

(16) 

C - MidP2 Intercept - LoMidP2DV - ( (UpMidP2DV - LoMidP2DV) 
/ (UpMidP2ML - LoMidP2ML) ) * LoMldP2ML 
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Substitution of these coefficients in formula (12), 
and multiplication by the scale factor to convert from pixels 
to millimeters, provides: 

5 (17) 

MLb - PixelSize * (MidP2Slope * Pt2ML - Pt2DV ♦ 
MidP2Intercept) / SQRT{MidP2Slope 2 ♦ 1> 

The same derivation provides the formula for the DV 
10 position of the second target point, defined as its distance 
within its section from the line between the earbars: 

(18) 

DV B - PixelSize * (BarbarSlope * Pt2ML - Pt2DV + 
15 * Barbar Intercept) / SQRT (BarbarSlope 2 + 1) 

The same reasoning leads to equations of the form 
of Equations (11), (17), and (18) to provide the positions of 
the first and third target points (ML^, DV A , AP A ) and (MI^, 

20 DV C , AP C ) in terms of millimeters from earbar zero. Note that 
no "BBZeroML" or "BBZeroDV" terms appear in (17) and (18) to 
match the "EBZeroAP" term in (11). This is because AP is 
defined as the position of a plane in a series of parallel 
planes, with a "starting plane 9 * based on the earbar entries, 

25 while ML and DV are defined specifically as distances from a 
point to a line in the same plane; for DV, of course, this 
makes the assumption that the earbar line would be in the same 
position in the three target-point slices as in the earbar 
slices. 

30 

6. A sixth step is simply to pass on the GDI- 
entered points (as converted to millimeter units and 
referenced to earbar zero) to the Core Algorithms outlined in 
Parts II and III. Other input options, however, are also 
35 useful. The most obvious of these is to enter coordinates at 
the starting point of th Core Algorithms, without reference 
to imaging data. This would be th approach, for example, if 
the user w re relying on a brain atlas, and this is expect d 
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to be a common approach for experimental applications 
involving rodents, which conform closely enough to a standard 
atlas, because of inbreeding, that imaging data on individuals 
may not be necessary. Other additional input options are 
5 described in Part VIII, "Alternative input methods" . 

One other variable to be entered by the user along 
with the original target points is SideArm, which indicates 
that the cannula-carrier is to be mounted on the right (-1) or 
left (+1) stereotactic arm. This variable of course has no 
10 effect on the Core Algorithms (Part II and Part III) , but it 
is important to the positioning of stereotactic equipment 
(Part V and Part VI) used as an example herein. 

7. A seventh step is to enter the physical 

15 characteristics of the coaxial elements (inner and outer 
diameters, cannula length) and collar mount elements 
(thickness) of the implantation tool. These parameters will 
not affect the trajectory itself, but they will affect the 
stereotactic positioning necessary to achieve the trajectory. 

20 The three original target points (Part I, Step 4) 

and the two tangent points (Part III, Step 1) on the circle 
are necessary and sufficient for calculating the trajectory. 
Additional details must be entered, however, to enable the 
accurate expression of the trajectory by the implantation 

25 tool, and to enable construction of the tool and the implanted 
compound tubing. The approach of the guide-tube cannula to 
the arcuate portion of the trajectory depends on accounting 
for a fine detail in the alignment of the curved tubing within 
the straight cannula. This detail is the fact that the 

30 curvature of the tubing (e.g., nitinol tubing) will begin to 
express itself slightly before the tubing emerges from the 
cannula. Mathematically, the only way to avoid this 
circumstance is to assume that the inner diameter of the 
cannula and the outer diameter of the nitinol tubing are 

35 identical, which is not physically possible. Since the 

tubing, therefore, must be smaller than the cannula, it is 
important to account f r the curvatur of the tubing that is 
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cone aled in th tip of the cannula. A comparable situation 
exists, of course, for the straight wire (e.g., nitinol wire) 
as it emerges from the curved nitinol tubing. These details, 
although small, can have a substantial effect on the actual 
5 trajectory expressed by the implantation tool. The following 
values are provided: 

LB - length of the guide- tube cannula 
10 * Inner diameter of the guide -tube cannula 
OD - outer diameter of the curved nitinol tube 
IDNitTube - inner diameter of the curved nitinol tube 
ODNitWire - outer diameter of the innermost nitinol wire 



10 



15 



(1) 
(2) 
(3) 
(4) 
(5) 



Finally, the thickness of the collar mounts and the 
20 thickness of the cannula-carrier block in which the guide-tube 
cannula (and therefore the implantation tool as a whole) is 
mounted are necessary for calculating the required lengths of 
tubing in the implantation tool: 
(6) 

25 CollarMountThickness - thickness of the collar mounts 

(7) 

Mount ingBlocfcThickness - thickness of the cannula- carrier block 

8. An eighth step is to make fine adjustments for 
30 the three initial target points. 

In practical use, for various reasons, it is 
helpful to be able to add slight adjustments to the initial 
target points, whatever input method is employed. These 
adjustments may be made at the outset or, more commonly, when 
35 the trajectory is being visually evaluated and readjusted 
(Part IV Step 3). 

FitRad, an optional small quantity to be added to 
th DVs of two targ t points, is a means of forcing the radius 



WO 96/03087 



PCT/US95/09391' 



28 

of the circle to a specific value. This conveniently permits 
matching the planned trajectory to the implantation tool 
available, for example. FitRad is an arbitrarily and 
empirically selected number that is added to the DV of target 
5 point 1 (DV A ) and to the DV of target point 3 (DV C ) . 

AvoidDivZero can be used to avoid a division-by- 
zero error for a trajectory that lies fully within a single 
coronal plane. AvoidDivZero is an arbitrarily and empirically 
selected very small number that is added to the AP of target 
10 point 1 (AP A ) and subtracted from the AP of target point 3 
(AP C ) . It can easily be made so small it has no practical 
effect on placement, but does permit the algorithms to proceed 
to calculate a very close approximation of the desired 
trajectory. 

15 TiltFactor makes it possible to move the exit point 

slightly, to avoid a critical structure at the surface such as 
the sagittal sinus when the wire emerges. TiltFactor is an 
empirically selected small number added to the ML of target 
point 3 (MLc) . 

20 

II. Core Algorithms: From three points to a series of circle 
points. 

25 1. A first step is to derive coefficients of the 

formula for the plane that contains the three original (non- 
collinear) target points (MI*, DV A , AP A ) , (MLb, DV b , AP b ) , and 
(MLc, DV c AP C ) from Part I. Again, there are a number of 
ways for accomplishing this task. One useful way is 

30 illustrated below. 

The general formula for a plane is: 
(1) 

35 H*X+I*Y + J*Z + K- 0 

(See, Love, * laments of Analytic Geometry MacMillan, NY, p. 
136, 1940) 
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The coefficients for this general formula that 
specify the plane containing the circle (the "circle plane") 
can be calculated from the three given points in a number of 
ways. One efficient way is illustrated below: 

Since all three points satisfy the general formula, 
one can substitute in (1) to get: 



(2) 



(3) 



(4) 



H * ML. ♦ I * DV A ♦ J * AP» + K - 0 



H * MLg + I * DV B + J * AP B + K » 0 



H * VILq + I * DV C + J * AP C + K - 0 



Subtracting (2) from each of the other Equations (1) , (3) , and 
(4) gives: 



20 (5) 



H • {X - ML A ) + I * <Y - DV A ) + J * (2 - AP A ) . 0 



25 



30 



(6) 



(7) 



H * (MLjj - ML A ) + I * (DV B - DV A ) + J * (AP B - AP A ) - 0 



H * (ML C - ML A ) + I * (DV C - DV A > ♦ J * (AP C - AP A ) - 0 



This homogeneous system of three equations gives 
rise to a determinant which is an alternative formula for a plane: 
(8) 



35 



(X - ML A ) (Y - DV A ) (Z - AP A ) 
(MLq - ML A ) (DV B - DV A ) (AP B - AP A ) 
(ML C - ML A > (DV C - DV A ) (AP C - AP A ) 



40 Evaluation f this determinant can be used to 

deduce a s t of valu s for the unknowns H, I, and J of the 
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homogeneous syst m. This is done by evaluating th minor 
related to each term in the first row (see Taylor supra, pp. 
522, 536): 



(9) 



10 



15 



20 



(10) 



(11) 



H - 



(DV f 



B 



DV A ) (AP B 
DV A ) (AP C 



AP A ) 



(MLq - ML A ) (AP B - AP A ) 

(ML C - ML A ) (AP C - AP A ) 

(MLg - ML A ) (DV B - DV A ) 

(ML C - ML A ) (DV C - DV A ) 



The fourth coefficient that corresponds to these is 
derived by substituting H, I, and J in Equation (2) and 
25 solving for K: 

(12) 

R - (H * ML A + I * DV A + J * AP A ) 

30 This is not, and does not need to be, a unique set 

of coefficients. Any multiple of the derived coefficients H, 
I, J, and K will be describing exactly the same plane 
(Equation (1)). 

35 2. A second step is to find the center of the 

circle through these three original points. 

Again, there are a number of ways to do this. One 
efficient way is illustrated below. 

Since the center of any circle lies on the 
40 perpendicular bisectors of any triangle on it, manufacture the 
perpendicular bisector planes of two of the edges of the 
triangle formed by the three original points, then solve the 
three equations (on for the plane containing the thr 
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original points and two for the perpendicular bisector 
planes) . 

The general formula for the distance (call it D) 
between two points (x 2 , y 2 , z 2 ) and (x 1# y 1# z x ) is: 

5 

(1) 

D - ± SQRT[ (x x - x 2 ) 2 + {y x - y 2 > 2 + (Zj - z 2 ) 2 ] 

The general formulas for the direction cosines of a 
10 line between two points (x 2 , y 2 , z 2 ) and (x lr y x , z x ) are: 

(2) 

A - <x x - x 2 ) / D 

(3) 

15 B - iy x - y 2 ) / D 

(4) 

C - (»! - * 2 > / D 

(See, Love supra, p. 129, and Protter et al.. College Calculus 
with Analytic Geometry, Addison-Wesley, Reading, HA, p. 538, 
20 1964.) 

The general formula for a plane, again, is: 
(5) 

25 AX + BY + CZ + E»0 

(See, Love supra, p. 136) 

For a plane perpendicular to a line, the direction 
cosines of the line may be taken as the first three 
30 coefficients in the formula for the plane (Love supra, p. 

137). Therefore the plane perpendicularly bisecting the line 
between the original target points (MI^, DV A , AP A ) and (MI*, 
DV B , AP B ) is described by the formula: 

35 (6) 

[ (ML A - MI*) / D ] * X ♦ [ (DV A - DVg) / D ] * Y ♦ 
[ (AP A - AP B ) /D]*Z+B-0 




WOWS/03087 



Since the plane passes through the midpoint between 
the two points, this midpoint may be substituted in Equation 
(6) to give: 

(7) 

[ (ML A - MLg) / D ] * (ML A ♦ MI^) /2 + [ (DV A - DV fi ) / D 1 * (DV A + 
DV B )/2 + [ (AP A - AP fi ) / D ] * (AP A + AP B ) /2 + B « 0 



10 



15 



Solving for E: 
(8) 

B - - H * I (ML A 2 - MLb 2 ) + (DV A 2 - DV B 2 > + (AP A 2 - AP B 2 ) ] / D 

Multiplying both sides of Equation (6) by -D, 
substituting E from (8), and rearranging gives: 



(9) 



(MLg - ML A ) * X + (DV B - DV A ) * Y + (AP B - AP A ) * Z 
M * [ (MLb 2 - MI^ 2 ) «• (DV B 2 - DV A 2 ) + (AP B 2 - AP A 2 ) ] 



20 Similarly, for points (MI^, DV C , AP C ) and (MI*, DV b , AP b ) , 



(10) 

(MLq - MI^) * X + <DV C - DV B ) * Y + ' (AP C - AP B > * Z - 
H * [ (MLq 2 - MI* 2 ) + (DV C 2 - DV B 2 ) + (AP C 2 - AP B 2 ) ] 

25 

The formula for the plane specified by the three 
points (coefficients given in Part II Step 1 Equations (9)- 
(12) ), and the two formulas ((9) and (10)) for the 
perpendicular bisector planes between two pairs of these 
30 points, can be solved simultaneously to give the location of 
the center of the circle through the three points. Again, 
there are a number of ways to do this. One efficient way is 
illustrated below. 
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First, a matrix is designated consisting of the 
coefficients of th three plan s: 



10 



15 



20 



»1 



Second, the planes are sorted to assure that a row 
with a leading zero (if it is present) is at the bottom, and 
another row with a leading zero (if it is present) is second 
up from the bottom. 

Third, elementary row operations are performed 
(Hoffman et al. f Linear Algebra . Prentice-Hall, Englewood 
Cliffs, NY, p. 6, 1961), each of which leaves the matrix row- 
equivalent to the original matrix, resulting in a row-reduced 
echelon matrix (Hoffman, supra, p. 12) . The right-hand 
column of values in the solved matrix represents the desired 
simultaneous solution of the three equations, and therefore 
the center of the circle, designated herebelow as (xx, yy, 
zz) . 



25 3* A third step is to find the radius of the 

circle. This can be easily calculated as the distance from 
any one of the points to the center: 



(1) 

30 rad - SQRTUxx - ML A ) 2 + (yy - DVa )2 ♦ {zz ' ^a' 2 ' 

4. A fourth step is to use the coefficients of the 
formula for the plane (H, I, J, K) , the coordinates of the 
center of the circle (xx, yy, zz) , and the radius of the 
35 circle (rad) to determine the Y coordinate of each of the two 
points of intersection of the circle with an XY plane at an 
arbitrarily selected position along the Z axis. 

Again, there are a number of ways to do this. One 
useful way is illustrated below. 
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Th DV coordinate of the intersection of the circle 
with a coronal plane at AP arctra j can be calculated by the 
formulas below, where: 

5 AA - H / SQRT( H 2 + I 2 + J 2 ) 

BB - I / SQRT( H 2 + I 2 + J 2 ) 

CC - J / SQRT( H 2 ♦ I 2 + J 2 ) 

DD - K / SQRT( H 2 + I 2 ♦ J 2 ) 

10 The formula for the circle plane is: 
(1) 

AA*X + BB*Y + CC*Z+DD-0 

The formula for the circle is: 
15 (2) 

(X - xx) 2 + (Y - yy) 2 + (Z - zz) 2 = rad 2 

The formula for the coronal plane is: 
(3) 

20 Z . AP arctraj 

Substituting (3) in (1) produces: 
(4) 

AA*X+BB*Y+CC* AP arctra j + DD - 0 

25 

Substituting (3) in (2) produces: 
(5) 

(X - xx) 2 + (Y - yy) 2 + (AP arctraj - zz) 2 - rad 2 

30 Solving (4) for X produces: 
(6> 

X - {-DD - BB * Y - CC * AParctraj > I 

Substituting (6) in (5) produces: 
35 (7) 

(t (-DD - BB * Y - CC * AP arctra j) / AA } - xx) 2 ♦ 
(Y - yy) 2 + <AP arctraj - ") 2 « rad 2 



Resolving (7) in terms of Y produces 
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(8) 



♦ b * Y + c 



where . 

5 a - (BB 2 / AA 2 + 1) 

b - (- 2 * yy + 2 * BB * (xx + <CC * AP arctra; : + DD) / AA) / AA) 

c - (yy 2 ♦ <<-CC * AP arctraj - DD) / aa - xx) 2 + (AP arctraj - zz> 2 - rad 2 ) 

Solving (8) for Y by using the quadratic equation 
10 gives two possible answers, corresponding to the DV coordinate 
at which the circle passes through the coronal plane above and 
below: 



(9) 

15 above: DV arctraj - Y - <-b + SQRT<b 2 - 4 * a * c) ) / (2 * a) 
(10) 

below: DV^^aj - Y - <-b - SQRT(b 2 - 4 * a * c) ) /{2 * a) 

20 5* A fifth step is to repeat the same calculation 

for each of a series of planes parallel to the arbitrarily 
selected XY plane. This can be accomplished simply by 
repeating the formula many times over, once for each coronal 
plane, over the range of the AP limits of the circle, making 

25 the calculations twice for each plane to account for the upper 
and lower intersection points. The selection of coronal 
planes may be made arbitrarily. It is convenient, however, to 
select the imaging scanner series of planes, at an interval of 
SliceThickness between selected planes, supplemented by the 

30 planes at the most posterior and most anterior points on the 
circle (see Part VII Step 4) . If an atlas is being used for 
entry of the original target points (see Part VIII Step 2) , 
then planes shown in the atlas should be included. 

6. A sixth step is to determine the X coordinate 

35 of each of the two points of intersection of the circle with 
each plane (in the series of XY planes parallel to the 
arbitrarily selected XY plan ) . Again, ther ar a numb r of 
ways t d this. One useful way is illustrated below. 
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Given that the AP of the coronal plan is 
designated AP arctraj , and the previously calculated DV is 
designated DV arctraj , the result follows simply from the 
substitution of these values in the formula for the circle 
5 plane: 



(1) AA*X + BB*Y + CC*Z + DD-0 

(2) 

10 ^arctraj - - CBB / AA) * DV arctraj - (CC / AA) * w ttCtraJ - (DD / AA) 

III. Core Algorithms: From a series of circle points and two 
tangent points to a series of trajectory points. 

15 

1. A first step is to select a point at which the 
straight guide-tube cannula will first contact the circle, and 
a second point at which the straight wire will first lose 
contact with the circle. These two tangent points must be 

20 entered by the user. 

There aire a number of ways in which the two points 
might be designated. One useful way is for the user simply to 
select the AP values of two of the parallel coronal planes for 
which the circle's intersection points have been calculated. 

25 These points are designated as (M^ tang , DV gtang , AP gtang ) and 
(Mlvtang^ ^wtang* AP^^g) , indicating their X, Y, and Z 
coordinates, respectively. 

2. A second step is to derive the direction 
cosines of the tangent to the circle, in the plane of the 

30 circle, at each of these two points. Again, the derivation 

could be performed in a number of different ways. One useful 

way is illustrated below. 

Since the tangent is perpendicular to the radius of 

the circle to its contact point, it must also lie in a plane 
35 normal to that radius. The formula for this plane is first 

d rived, th n the formula for the plane of the circl is 
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restat d, and the two formulas are solved simultaneously to 
give direction cosines that specify th tangent line. 

The radius of the circle that reaches the contact 
point of the tangent is a line, and therefore has direction 
5 cosines derivable from the coordinates of two points (x x , y 1# 
z x ) and (x 2 , y 2 , z 2 ) on the line and the distance d between 
them (Pr otter et al., supra, p. 538, and Love, supra, p. 129): 

(1) 

AlphaDi recCos - (x 2 - x x ) / d 

10 (2) 

BetaDirecCos - <y 2 - Y\) / d 

(3) 

GamnaDirecCos » (z 2 - z A ) / d 

15 Substituting the designations for the center of the 

circle (xx f yy, zz) , for the tangent point (MLg tang , DV gtang , 
AP gtang ) , and for the distance between then (SQRT((xx - 

W^tang) 2 + (YY ~ DVgtang) 2 + < 2Z " ^gtaog) 2 )) renders the 
equations to: 
20 (4) 

AlphaDirecCoBRftjjjus - (xx - MI^ tang ) / SQRT( (xx - MI^^) 2 ♦ 
(yy - DV g tang> 2 + <■« " AP gtang )2 > 

(5) 

25 BetaDirecCoSRADHjg - (yy - DV gtang ) / SQRT( (xx - MLg tang ) 2 ♦ 

(yy - DVgtang) 2 + (ZZ - APg tan g) 2 ) 

(6) 

GammaDirecCoB^^jug « (zz - AP gtang ) / SQRT( (xx - MLg tang ) + 
30 (yy - DV^^) 2 ♦ <»« " ^gtang^ ) 

The tangent lies in a plane normal to this radius. 
The general formula for a plane passing through a given point 
(x 0 , Yo» z o) and perpendicular to a given line (direction 
35 cosines A, B, C) is (Protter et al., supra, p. 545; and Love, 
supra, p. 136) : 
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(7) 

A * (X - x Q ) ♦ B * (Y - y 0 ) + C * <Z - 2 0 ) - 0 

5 Substituting in Equation (7) the point (111^.^, 

DV gtang' Mgtang) and the direction cosines of the radius given 
in Equations (4) , (5), and (6) provides an equation for the 
plane normal to the radius: 

10 (8) 

I {xx - MLg tang ) / SQRT( (xx - MLg tang > 2 + (YY - Wgtang )2 ♦ 
(zz - AP gtang ) 2 ) ) * (X - MI^ tang ) + t (yy - DV gtang ) / SQRT( (xx - 

MLgtang) 2 * <™ ' DV gtang> 2 + {zz ' ^gtang* 2 > 1 * < y " Dv gtang> ♦ 
[ (zz - AP gtang ) / SQRT((xx - Ml^tang) 2 + (yy - DV^^g) 2 + (zz -AP gtang > 2 

15 ) ] * (Z - AP gtang ) - 0 

Nov, the tangent line lies not only in the plane 
normal to the radius (Equation (8)), but also, by definition, 
in the circle plane: 

20 

(9) 

AA*X + BB*Y + CC*Z +DD- 0 

All solutions to the two Equations (8) and (9) have 
25 the proportions of the determinants implicit in the equations 
(Taylor, supra, pp. 522, 546, 554). For illustration, the 
equations and the determinants can be restated in more general 
form: 

30 (10) 

A*X + B*Y + C*Z + D«0 

(ID 

A' * X ♦ B' * Y + C * Z ♦ D' - 0 



35 
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The proportional determinants implied by the pair 
of equations are: 
(12) 



B C 




C A 




A B 


B'C 




C f A' 




A'B' 



10 



IS 



20 



25 



30 



35 



The sane values can be expressed algebraically 
according to the standard definition of a 2x2 determinant: 
(13) 

(B * C - B' * C) : (C * A' - C * A) : (A * B' - A' * B) 

These are the proportions of the direction cosines 
of the intersection line that is contained in both planes (10} 
and (11) . By substitution in (13) of the values of the 
coefficients of the plane of the tangent (normal to the radius 
to the tangent point) (Equation 8) and the plane of the circle 
(Equation 9) , and rearrangement, formulas for the direction 
cosines of the tangent can be expressed as three fractions, 
each of which has the same denominator: 

den - SQRT(P + Q + R) 

where 

P - <(<DV gtang - yy)) * CC - <<AP gtang - 22)) * BB)* 

Q - ( ( (AP^^ - M>) * AA - ( (ML^ tan g - xx)) * CC) 

R - ( { (MLwf.. na - xx) ) * BB - ( (DV^. M „ - yy) ) * AA) 



(14) 
(15) 
(16) 



AlphaDirecCos gtang - ( (DVg^g - yy) * CC 
BetaDirecCosgtang - ( (AP, 



( « DV gtang " Vy>> * *A> 2 

gtang " "> * ra >/ den 



'gtang * «■> * " lM Vang " xx) * CC)/ den 
OamnaDirecCoSgtang - { (MLg tang - xx) * BB - {UV^^ - yy) * AA) / den 



40 



The actual value of the square root in the 
denominator can be either positive or negative. AA, BB, CC, 
and DD (as d fined earlier) represent th direction cosines 
for th circle-plane — i.e., for the vector normal to the 
plane. The dir ction of this vector depends in part on the 
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10 



ordering of the coronal planes (Z = AP A , Z = AP B , Z = AP C ) 
containing the three original targ t points. When these 
planes proceed from posterior to anterior — from smaller to 
larger AP — the normal is directed to the "subject's right" 
half of the coordinate space (octants II, III, VI, VII); when 
the ordering is from anterior to posterior — from larger to 
smaller AP— the normal is directed to the "subject's left" 
half of the coordinate space (octants I, IV, V, VIII) . The 
direction does not affect the calculation of the circle's 
intersection points with coronal planes, based on AA, BB, CC, 
and DD (see above) , but it does affect the direction cosines 
of the guide-tube cannula and the wire, which in turn affect 
the orientation of the guide-tube cannula. A change in the 
sign of "den", changing the sign of the direction cosines, 
15 reverses the direction of the cannula. When the order of AP A , 
AP B , AP C proceeds from anterior to posterior (AP A >AP C ) , a 
positive den directs the guide-tube cannula from its own top 
to its own tip, so that if the tangent point is posterior to 
the lowest point on the circle, then DVtop is positive and 
20 DVtop-DVtip as calculated below (Part V Step 2 Equation (8)) 
is positive; when the order of AP A , AP B , AP C proceeds from 
posterior to anterior (AP C >AP A ) , a negative den has the same 
effect. Therefore, in order to permit the ordered input of 
points AP A , AP B , and AP C in sequence either from anterior to 
25 posterior or from posterior to anterior, SIGN(AP A - AP C ) is 
factored into den. This has the effect of directing the 
guide-tube from its own top to its own tip in all cases. If 
the tangent point is posterior to the lowest point on the 
circle, then DVtop is above DVtip and DVtop-DVtip is positive; 
if the tangent point is anterior to the lowest point on the 
circle, then DVtop is below DVtip and DVtop-DVtip is negative. 
The complete formula for the denominator is therefore: 
(17) 

den - <SIGN(AP A - AP C )) * SQRT( 



30 



35 <nDV gtang -yy)) * CC - ( <AP gtang -M> ) * BB) ♦ 

l<<AP gtang -">> * ** ■ <«%tang-**» * * 
(((M4 tang -xx)) * BB - <<DV gtang -yy>> Afc) 2 ) 
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An entir ly analogous derivation leads to th 
dir ction cosines of the wire tangent AlphaDirecCos^ang , 
BetaDirecCos^ang, and GammaDirecCos^^g. 

3. A third step is to use the direction cosines to 
5 determine the points at which the guide-tube cannula and the 
wire intersect each plane (in the series of planes parallel to 
the arbitrarily selected plane) . Again, there are a number of 
ways to do this. One useful way is illustrated below. 

Designate the coordinates of the point of 
10 intersection of the tangent with the parallel plane as 

M^egtraj X-coordinate) , DV legtraj (the Y-coordinate) , and 

APiegtraj I th * Z-coordinate, which is the Z-coordinate of the 
entire parallel XY plane (coronal plane) in this example) . 

Mow the direction cosines, of a line, as stated 
15 above, are related to the coordinates of any two points (x t , 

y lf z x ) and (x 2 , y 2 # * 2 ) through which the line passes, and the 
distance d between the two points (Protter et al., supra, p. 
538, and Love, supra, p. 129): 
(1) 

20 AlphaDirecCos - (x 2 - »i) / d 

(2) 

BetaDirecCos - <y 2 - y^) / d 

(3) 

GamraaDirecCos - (» 2 - *i> / d 

25 

It is evident, from dividing (1) by (3), that 
(4) 

AlphaDirecCos / GammaDirecCos - (x 2 - x x ) / (a 2 - 

30 The line passes -through the tangent point (KLg tang , 

DV^.^, APnt-Ma) and through the point in the parallel plane 

(W^egtraj* DVlegtraj' A*legtraj> . °* ^ icn ^legtraj <«»• 
Z-coordinate) is known. Therefore, by substitution of these 
points for the two points (x 1# y lr z x ) and (x 2 , y 2 , x 2 ) in (4), 
35 it is clear that: 

AlphaDirecCos gtang / GanrnaDirecCoB gtang - <ML legtra j - M^gtang 5 ' 
(AP legtraj " APgtang* 
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and by r arrangement, 
(6) 

^legtraj " ^tang + (AlphaDirecCos gtang / GammaDirecCos gtang ) * 
^ AP legtraj " ^gtang* 

5 

In similar fashion, it is clear from this that: 
(7) 

DV legtraj - ^tang + <BetaDirecCos gtang / GammaDirecCos gtang ) * 
(AP legtraj " AP gtang ) 

10 

The guide-tube cannula tangent coordinates (Ml^ tang r 
DV gtang , APg t ang) can of c° urse replaced by the analogous 
wire tangent coordinates (KL^^g, DV^^g, AP^^g) , in which 
case the resulting intersection coordinates (ML legtra j, 
15 DV legtraj , AP legtraj ) describe the wire tangent trajectory. 

4. A fourth step is to construct a series of X, Y, 
and Z values copied from the results of Part II Steps 4-6 and 
Part III Step 3, which provide the points on the curved and 

20 straight parts of the trajectory, respectively. On the basis 
of the two tangent points selected in Part III Step 1, this 
series contains the coordinates of the guide-tube cannula 
intersections with the parallel planes up to the point of 
tangency, then the coordinates of the circle intersections up 

25 to the second point of tangency, then the coordinates of the 
wire intersections for the remainder of the parallel planes. 
This series represents the required trajectory, QED. 



30 IV. Machine-dependent output Algorithms: Prom a series of 
trajectory points to a visual display. 

1. A first step is to convert the final output 
trajectory (Core Algorithms) into pixel positions in a TIF 
35 file to be displayed by an imaging system. There are various 
ways to accomplish this conversion. One reasonable example is 
given b low. 
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The AP c rdinates for a trajectory point, 
(AParctraj) or (^legtraj) » provide th equivalents in TIF 
coordinates (AParctj-ajpuj) or (AP legtrajpIX ) directly from Part I 
Step 5 Equation (11) . 

5 

(1) 

AP B m SliceThickness * (Pt2AP - EBZeroAP) 

Substituting AP arctraj for AP B and AP arctrajPIX for Pt2AP gives: 
10 (2) 

AP arctraj " SliceThicJcness * (AP arct ra j p ^ - EBZeroAP) 

and rearranging provides: 
(3) 

15 ^arctrajPDC " EBZeroAP ♦ AP arctraj / SliceThickness 

The ML and DV coordinates for a trajectory point 

(•tt-arctraj anA ^arctraj' or ' (^egtraj and ^legtraj)' 1x1 terms 
of the orthogonal coordinate system used by the Core 
20 Algorithms, can be scaled and converted to pixels in the TIP 

coordinate system. Equations (17) and (18) from Part I Step 5 
are already available for relating MRI ML and DV values in 
pixels (such as Pt2ML and Pt2DV) to Core Algorithm values in 
millimeters (such as MLq and DV B ) : 

25 

(4) 

MLq - PixelSize * (MidP2 Slope * Pt2ML - Pt2DV + MidP2 Intercept) / 
SQRT(MidP2Slope 2 + 1) 

30 (5) 

DV B - PixelSize * (BarbarSlope * Pt2ML - Pt2DV + Barbar Intercept) / 
SQRT {BarbarSlope 2 + 1) 

The corresponding equations for trajectory points 
35 on the circle can be written by substituting the known values 
ML arctraj and DV arctraj ( Part II* Step 4, in millimeters) for 
MLg and DV B respectively in (1) and (2) , and the unknown 
values ML^^^jpjx and DV arctra j PIX for Pt2ML and Pt2DV: 
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(6) 

^arctraj " PixelSize * (MidP2Slope * ML arctrajpix - DV arctrajpIX + 
MidP2 intercept) / SQRT(MidP2Slope 2 + 1) 

5 

(7) 

^arctraj - PixelSize * (EarbarSlope * ML arctrajpIX - DV arctrajpIX + 

Barbar Intercept) / SQRT( EarbarSlope 2 + l) 

10 The only unknowns in these two equations are 

ML arctrajPIX and DV^^p^, Points in MRI pixel terns that 

are needed. Solving the two equations for the two unknowns 
provides the formulas: 

15 (8) 

^arctrajPIX " 

[ (DV arctra j / PixelSize) * SQRT (EarbarSlope 2 ♦ 1) 

- (ML arctraj / PixelSize) * SQRT ( Mi dP2 Slope 2 + 1) 

+ MidP2 intercept 

20 - Barbar intercept ] 

/ [ EarbarSlope - MidP2Slope 1 



(9) 

^arctrajPIX " 

25 I (DV^^aj / PixelSize) * MidP2Slope * SQRT (EarbarSlope 2 + 1) 
_ (ML^traj / PixelSize) * EarbarSlope * SQRT(MidP2Slope 2 ♦ 1) 
+ EarbarSlope * MidP2 Intercept - MidP2Slope * Barbar Intercept ] 
/ I EarbarSlope - MidP2Slope ) 

30 The same derivation provides the analogous display- 

ready leg-points (ML legtra j, DV le gtraj' A p legtraj) • 

2. A second step is to perform operations similar 
to those of Part I Step 2, configuring equipment that can be 
used to interface between the analytic algorithms that plan 

35 the trajectory and the display from the imaging scann r. 
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Again, there are a number of ways to do this. On useful way 
is illustrated below. 

As explained above (Part I Step 2 Procedures (1) to 
(3)), the TIF file as converted for video display is a series 
5 of numbers (of magnitude between 0 and 254) representing gray- 
scale values of pixels in an orderly array. For displaying 
point (ML arctrajPIX , DV arctrajPIX , AParctrajPix) • Procedures (5) 
and (4) of Part I Step 2 are reversed: 

(1) The analytic algorithms pass this converted 
10 trajectory point to a text file, which is then taken up by 

another 80486 program (in Microsoft Visual BASIC) that surveys 
the original TIF file from the MRI and highlights the pixel 
representing the trajectory point in the appropriate section 
(AParctrajPix) • 11118 Program uses the TIF file's row and 

15 column arrangement and its header length to determine which 
pixel position to highlight. Altogether there are 
(HeaderLength + 256*256) pixels. Actually pixels in a 3x3 
square grid are highlighted (set to 255) for each point. A 
slightly different shape is used for the trajectory points 

20 nearest the original target points, and other shapes are used 
for the points at the beginning and end of the "target 
structure" and for the entry and exit points on the cranium 
(see Part VII Step 6 below) • 

(2) The file is then displayed by a commercial 
25 image-display package (Adobe Photoshop) , for visual 

observation. The highlighted pixels appear distinctly, at the 
predicted anatomical position of the point. 

3. EPRON: A third step is to evaluate the 
30 trajectory resulting from the entry of the three original 

target points and two tangent points. If visual inspection of 
the display shows that the anatomical placement could be 
improved, changes can be made. The numerical value of an 
original target point or points can be modified. The 
35 algorithms repeat the analysis and display a modified 

trajectory. A sophisticated system (St ps 3b and 3c below) 
may be implemented utilizing procedur s and s ftware stored on 
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a perman nt memory device , such as an erasabl programmable 
read only memory (EPROM) chip programmed to make these 
modifications by a variety of GUI input options and to keep 
all other algorithmic computations updated accordingly. 

5 

a. Real Input Mode 

An elementary way to modify the input is simply to 
switch to "Real Input 9 * mode (see Part I, Step 6) , enter 
modified target points, and allow the algorithms to repeat the 
10 analysis. The specific "fine adjustments 1 * introduced above 
(Part I Step 8) may also be made at this point. 



b. Cross-section manipulation 

Another way will be to incorporate image display 

15 algorithms into an EPROM device to permit rapid processing of 
GUI input and output. Both the underlying calculations and 
the image displayed on the basis of the calculations will be 
rapidly updated in response to keyboard or GUI inputs. The 
EPROM will display a series of separated parallel cross- 

20 sections that can be manipulated in perspective. The cross- 
sections can be viewed from an angle, with the dorsal edge up 
and the ventral edge down (for brain sections) and with the 
anterior surface appearing in perspective. Each trajectory 
point appears in position in its section. By option, the 

25 display can be adjusted, and the cross-sections can be 

manipulated in several ways by moving a cursor, e.g., with a 
mouse: 

(1) the cross-sections can be variably rotated, 
each on its own mid-sagittal axis, like slats of a blind? 
30 (2) the screen distance between cross-sections in 

the display can be changed; 

(3) cross-sections can be added to or removed from 
the display; 

(4) a common portion of each cross-section can be 
35 erased from the display for better visibility; 
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(5) the row of cross-sections can be divided at any 
point leaving a single larger between-section interval between 
two groups of cross-sections; 

(6) a given cross-section can be removed from the 
5 row and examined full-front; 

(7) that cross-section or all the cross-sections 
can be magnified or diminished; 

(8) alternatively, the screen can be split between 
a smaller area for a reduced rendition of the series of cross- 

10 sections, emphasizing the trajectory points, and a large area 
showing selected cross-sections full -front. 



c. Trajectory manipulation 

The display of trajectory points can be separately 

15 manipulated. Its portrayal on the screen can be shifted from 
single points on individual sections to a continuous line of 
points for manipulation as a unit under optional constraints. 
The basic pattern of interaction is simple. The user "clicks 
and drags" a single point on the line from one position to 

20 another, the line remodels itself in response, and all 

displays and algorithms are updated. A given version can be 
saved at any point, to be recalled later on. Some of the 
constraints are: 

(1) the entire trajectory moves as a single object, 
25 without changing its curvature or its angles of orientation; 

(2) the line is fixed at two selected (clicked) 
points and when a third point is dragged it rotates the 
trajectory around these points as an axis; 

(3) either the guide-tube cannula tangent or the 
30 wire tangent is selected as an axis, and when a point is 

dragged the trajectory rotates on that axis, without otherwise 
changing; 

(4) each tangent is fixed at one selected point, 
and when a point on the trajectory is dragged the curvature of 

35 the arc is increased or decreased, and tangents are realigned 
to remain tangential, while the angles f orientation of the 
circle plane remain the same; 
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(5) the angles (direction cosines) of the tangents 
are fixed, and when a point on the arc is dragged the 
curvature of the arc is increased or decreased, and tangents 
are translated to remain tangential, while the angles of 
orientation of the circle plane remain the same; 

(6) the point selected is substituted for the 
nearest original target point or tangent point, and the 
trajectory is adjusted in response, without altering the other 
four entry points. 



A digital display of critical features of the 
trajectory will also be available, and changes made on this 
display from the keyboard will be reflected in the display of 
the trajectory. These features will include the coordinates 

15 of the original target points and tangent points, the center 
and radius of the circle and the angles of orientation of its 
plane (or the various features of the helix listed in Part IX, 
Step 15) , the angles of orientation of the tangent lines, the 
coordinates of the points at which the trajectory enters or 

20 exits the body, and the coordinates of the points at which the 
trajectory enters or exits the target structure. 

Input options will also include the helical 
trajectory, and the complex trajectories consisting of 
alternating straight and curved portions (see Fig. 2A-2E) . 

25 More than three target points can be entered for fitting to 

these trajectories. The user will be required to specify the 
number of full periods to be included, for a helical 
trajectory^ and the locations of break-points between one 
curvature and another for a complex trajectory. Options 

30 comparable to those listed above will be included for these 

trajectories in the EPROM. The arcuate trajectory will still 
be available by option, and the user will be able to shift 
back and forth among the best-fitting arcuate, helical, and 
complex trajectories for the full set of target points 

35 entered. 



WO 96/03087 



PCT/US95/09391* 



49 

V. Machin -dependent: output Algorithms: Prom th planned 
position of the guide-tube cannula to prescribed stere tactic 
angular settings. 

5 1. The first step is to calculate lengths 

associated with the curvature concealed in the tip of the 
guide-tube cannula 20 and the straight wire segment 24 
concealed in the tip of the curved nitinol tubing 28. There 
are a number of ways to do this. The example below describes 

10 one technique, referring to Pig. 5, where a partial cut-away 
view of the tool of Fig. 1 is shown. 

The length of the guide-tube cannula 20 in which 
resides the initial curvature of the curved nitinol tubing 28 
can be calculated by considering that this curvature is 

15 restrained by the inner edge of the cannula 20 at the cannula 
tip (point P) . It is assumed that the curved nitinol tubing 
28 is straight within the lower part of the cannula, held 
against its wall, just above a point at a level N near the 
tip, at which point it curves down and across the dead space 

20 in the cannula, from the far wall (far from the crater of the 
circle) to the near wall (near the center of the circle) . It 
emerges pressed against the near wall at the tip (point P) . 

The line from the center (point C) of the circle to 
a point £ on the cannula directly across the cannula from 

25 point N is assumed to be perpendicular, since the tubing is 
assumed to express its curvature as far up the cannula as 
possible. The length of such a perpendicular to the inner 
edge of the curved tubing (point N) is simply the radius (rad) 
of the circle. (At this point the precise definition of the 

30 radius as the distance from the center to the inner aspect of 
the curved tubing becomes geometrically important and 
desirable from a practical point of view, in that it is easily 
measured and related in the simplest way to the fabrication of 
the curved tubing by wrapping it around a mandrill • ) Its 

35 length to E is rad, less the deadspace (the difference between 
the inner diameter of the cannula and the uter diam ter of 
the curved tubing) . The distance from th center of the 
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circle to the curved tubing at its point of emergence P, where 
it is pressed against the inner wall of cannula, is rad. 

In aCEP, the leg of the triangle EP, along the 
inner edge of the cannula, is the length (CurvelnGuideLength) 
5 of the terminal portion of the cannula within which the 
initial curvature of the circle is contained. It can be 
calculated from the Pythagorean theorem (CE 2 + EP 2 = PC 2 ) 
that: 

10 (1) 

(rad - (ID - OD)) 2 + CurvelnGuideLength 2 * rad 2 

or, 
(2) 

15 CurvelnGuideLength - SQRT(rad 2 - (rad - (ID - OD)) 2 ) 

On the same triangle, the angle ECP that subtends the arc 
concealed in the cannula can be calculated as: 

20 (3) 

arcsin (CurvelnGuideLength / rad) 
and therefore the linear length along this arc must be: 
25 (4) 

CurvelnGuideArcLength • rad * arcsin (CurvelnGuideLength / rad) 

At the tip of the curved nitinol tubing, the forced 
curvature of the natively straight nitinol wire 24 will end at 

30 a point T on the inner wall of the tubing slightly before its 
tip. In the final arc ("omega") of the tubing, the wire is 
straight, traversing the tubing deadspace. At the tip of the 
tubing, the radius from the center of the circle (C) to its 
inner wall (H) does not quite reach the edge of the wire. It 

35 has to be extended across deadspace. But this deadspace is 
not simply the diff rence between th diameters, because the 
wire crosses the lumen at the tip of the tubing at an angle. 
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By similar triangles (aTCX and aXFR) , however, it can be seen 
that the angle XFR is also omega. Th refore: 

(5) 

5 FX - RF / cos (omega) - ODNitWire / cos (omega) 

Since HF = IDNitTube, the length HX can now be written: 

(6) 

10 HX - IDNitTube - (ODNitWire / cos (omega)) 

and CX, the hypotenuse of aTCX, is (rad + HX) , or: 

15 rad + IDNitTube - (ODNitWire / cos (omega)) 

The cosine of the angle omega in aTCX can now be set equal to 
CT/CX, leading to: 

20 (8) 

cos (onega) - rad / (rad + IDNitTube - (ODNitWire / cos (omega)) 
Solving (8) for cos (omega) yields: 
25 (9) 

cos (omega) - (rad + ODNitWire) / (rad + IDNitTube) 
Therefore omega can be identified: 
30 (10) 

omega - arccos ((rad + ODNitWire) / (rad + IDNitTube)) 

and the linear length of the arc subtended by omega, the 
length, that is, of the tubing that conceals straight wire, 
35 is: 

(11) 

StraightWireArcLength - rad * arccos ((rad ♦ ODNitWire) / (rad + 
IDNitTube) ) 



■ . . . 
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7^ 



The length (SF) of the straight wire "concealed" in the tubing 
(although aXFR is beyond the tip) is also easily calculated 
from omega/ since: 



5 (12) 

tan (omega) - SF / CS 
or, substituting and rearranging, 



10 (13) 

SP - (rad + ODNitWire) * tan (omega) 
and, substituting the original values: 
15 (14) 

StraightWireLength - (rad + ODNitWire) * tan( arccos ((rad + 
ODNitWire) / (rad + IDNitTube) ) ) 

The length (GTUpper) of the guide-tube cannula 20 
20 down to the beginning of the curvature of the nitinol tubing 

is a significant value, because it represents the total length 
of the guide-tube cannula line to its tangent point, and is 
used in certain calculations below: 



25 (15) 

GTUpper - I*E - CurvelnGuideLength 

2. A second step is to calculate the locations of 
the top and the tip of the guide-tube cannula 20 when it is in 
30 its final position. Again, there are a number of ways to do 
this. One useful way is illustrated below. 

As stated above (Part III Step 2, Equations (1-3)), 
the general formulas for the direction cosines of a line 
between two points (x 2 , y 2 , z 2 ) and (x lf y x , z x ) at a distance 
35 d are: 
(1) 

AlphaDirecCos « (x 2 - x x ) / d 

(2) 
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BetaDirecCos - (Y2 - Xi> / d 

(3) 

GammaDirecCos - <z 2 - z x ) / d 

(Love, supra, p. 129, and Protter et al., supra, p. 538) 

5 

The point of tangency for the guide-tube cannula 
(MVang* w gtang* AP gtang ) is known and can be substituted for 
one point, and the coordinates for the top of the guide-tube 
cannula (MLtop, DVtop, APtop) can be substituted for the other 
10 point in these formulas. The distance d, furthermore, from 
the top of the guide-tube to the point of tangency (at which 
curvature begins) , is now known as (GTOpper « LE - 
CurvelnGuideLength) , derived above. 



15 (4) 
(5) 



AlphaDirecCos gtang - (MLgtang ' MLtop) / GTOpper 

BetaDirecCos gtang - ^gtang " Wtop) / GTOpper 

(6) 

20 GainmaDirecCoSgtang - (AP gtang " APt °P ) / GTOpper 

From these formulas, by rearrangement, it 
immediately follows that: 





(7) 










25 


(8) 


MLtop 


" ^^gtang 


- GTOpper 


* AlpbaDirecCo8 gtang 






DVtop 


" ^gtang 


- GTOpper 


* BetaDirecCoSgt^g 




(9) 














APtop 


" ^gtang 


- GTOpper 


* Gan*naDirecCos gtang 



30 

In like manner, the coordinates for the tip of the 
guide-tube cannula (MLtip, DVtip, APtip) can be substituted in 
formulas (1-3), along with the coordinates of the point of 
tangency and its distance from the tip (CurvelnGuideLength) , 
35 to provide the coordinates of the tip of the guide-tube 
cannula: 
(10) 

MLtip - MLgtang - CurvelnGuideLength * AlphaDirecCos gtang 
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DVtip 

(12) 

APtip 

5 

3. A third step is to generate the angles at which 
to set the stereotactic apparatus. There are a number of ways 
to do this, depending in part on the type of stereotactic 
apparatus to be employed in surgery. The example below 

10 describes one technique adapted to the Kopf stereotactic 

apparatus. All angles have been left in radians, although in 
the software itself they must be converted into degrees to 
conform with the units used in typical stereotactic scales. 

One specific implementation of an implantation tool 

15 10 according to the present invention is shown in Fig* 4, 

where a Kopf stereotactic device 100 is shown with a cranium 
150 mounted thereon. Briefly, the device 100 is used to 
position the tool 10 by holding the cranium 150 in a fixed 
position and by moving the guide-tube cannula relative to the 

20 cranium 150 through one of several calibrated adjustments. 
The device 100 includes horizontal and vertical slides 195, 
120 and two adjustment scales 160, 170. The scales are known 
as the PED scale 160 and the ROT scale 170. The cranium 150 
is fixed through use of an earbar 180 and eyebars 140. The 

25 implantation tool 10 is fixed to the device 100 by inserting 
the guide-tube cannula 20 into a block 15 mounted on a rod 
held by the horizontal slide 195. Calibration of the device 
100 for use with the tool 10 of the present invention will now 
be described. 

30 The Kopf stereotactic apparatus provides for two 

calibrated rotational adjustments. The PED scale 160 is 
located at the joint between the cannula-carrier and the 
lateral arm 110 of the apparatus on which it slides, and 
permits a rotation of the carrier around the DV axis by a 

35 measured angular distance (PED) away from the sagittal plane 

(YZ plane) . Th physical axis f rotation for this adjustment 
is ffset fr m the shaft of th guide-tube cannula 20 itself, 
in the Kopf st reotactic, but this offset can be compensated 



° DV gtang " Curve InGuideLength * BetaDirecCoSg tang 
- APg tang - CurvelnGuldeLength * GammaDirecCoSg tan g 



WO 96/03087 



PCI7US95/09391 



55 

by a translational adjustment lat r on. In effect, adjustment 
at this joint rotates what will be called the "cannula-carrier 
plane 1 * away from its initial orientation parallel to the 
sagittal plane (YZ plane) of the cranium 150, to a new 
5 orientation, parallel to the plane containing the guide-tube 
cannula 20 in its planned final position and the projection of 
the guide-tube cannula onto the axial plane. This new 
position is specified by the direction cosines of the guide- 
tube cannula, which imply a projection onto the axial plane 

10 (XZ plane) at an angle PED from the Z axis. 

The angle FED can be calculated from the known 
points at the top and the tip of the guide-tube cannula by 
considering this projection (Pig. 6A) . It is a transverse 
line in the axial plane that constitutes the hypotenuse of a 

15 triangle of which the height (parallel to the AP axis) is: 

(1) 

A - APtop - APtip 

20 and the base (parallel to the ML axis) is: 
(2) 

B - MLtop - MLtip 

25 The length of the hypotenuse is therefore: 
(3) 

C - SQRT( (APtop - APtip) 2 + (MLtop - MLtip) 2 ) 

30 The angle of the hypotenuse from the AP axis is therefore: 
(4) 

PED unadju8ted ~ arccos < (APtop - APtip) / SQRT( (APtop - APtip) 2 + 
(MLtop - MLtip) 2 ) ) 

This formula for the angle PED needs to be 
adjusted, howev r, to accommodate the ambiguity f the arccos 
functi n and th peculiariti s of the stereotactic equipment. 
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For each cosine value (from 1 to -1) the spreadsheet arccos 
function choos s the corresponding arccos valu betw en 0 and 
7r although there is also a distinctive angle between rr and 2n 
which should be chosen in some cases. Fig. 6A with labeled 
5 vertices illustrates this point, portraying the projection of 
the guide-tube cannula 20 onto the XZ plane, translated to 
bring its tip to the origin for convenient analysis. The 
angle PED^^jugt^ specified by Equation (4) is the angle ROG, 
measured from the positive AP axis to the projection of the 

10 guide-tube cannula. The arccos function returns only the 
arccos value between 0 and tt, even though in the case 
illustrated the angle ROG would be properly designated by the 
corresponding angle between ir and 2 it, rotating through 
quadrants I (+AP, +ML) and II (-AP, +ML) and into Quadrant III 

15 (-AP, -ML) . One must account for the full 360° range of 

possible angles, counting in order from the "Anterior" axis 
through quadrants I, II, III, and IV. As the cosines go from 
1 to -1 through quadrants I and II — i.e., (MLtop - MLtip) > 0 
— the angles provided by the spreadsheet go from 0° to 180°. 

20 As the cosines continue from -1 to 1 through quadrants III and 
IV, however — i.e., (MLtop - MLtip) < 0 — the angles provided 
by the spreadsheet do not continue from 180° to 360° as 
desired, but go from 180° back down to 0°. This second series 
must be corrected by taking 360 - PED unadjuBted' instead of 

25 PED unadjusted» 

The angle PBDunadjusted of Equation (4) also needs to 
be adjusted to accommodate peculiarities of the stereotactic 
equipment. The Kopf stereotactic scale for PED has a 
conveniently readable excursion of only ±50°. In keeping with 

30 the convention that anterior is positive, it is consistent to 
express the adjustment as a positive angle when it shifts the 
movable line on this scale toward the anterior, and as a 
negative angle when it shifts the movable line toward the 
posterior. This has to be true regardless of whether the 

35 carrier is placed on the left arm (SideArm -= 1) or the right 
arm (SideArm = -1) . 
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An analysis of the possible cases is necessary to 
generate the appropriat adjustment for a PED angl that 
rotates the cannula-carrier plane into each quadrant of the 
ML/AP plane. In Table 1, on the next page, ORI represents the 
5 original unadjusted angle PED unad j usted of Equation (4). 
Further simplification yields the final Table 2 on the 
subsequent page. (The variable names M APorient N , "MLorient" , 
and "SideArm" are introduced for SIGN (APtop - APtip) , SIGN 
(MLtop - MLtip) , and Carrier on the Right or Left, 

10 respectively. The "Total Adjustment* column is simplified, an 
example is introduced, and the row order is rearranged.) The 
eight combinations of sign conditions represented in these 
eight rows can be summarized in the formula that heads the 
right-hand column, because: (1) If APorient is negative, then 

15 180 must be added or subtracted. (2) If SideArm * NLorient is 
negative, the final value is either (ORI) or (ORI - 180) • 
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The appropriately adjusted formula is th r fore: 

(5) 

5 PHD - {- SideArm * MLorient) * ( [PED unad j usted ] - (180 * 

(-APorient)) ) 

or, 
(6) 

10 PED - (- SideArm * MLorient) * { tarccos ( (APtop - APtip) / 

SQRT( (APtop - APtip) 2 + (MLtop - MLtip) 2 )) ] - <180 * (-APorient)) ) 

Note that the set of cases in the table above 
assumes always an approach in which DVtop > DVtip. The 

15 alternative, an approach from below, is not permitted by the 
Kopf stereotactic equipment. 

The Kopf stereotactic apparatus also provides for a 
second calibrated rotational adjustment ("ROT scale") 170, 
permitting rotation of the guide-tube 20 within the originally 

20 parasagittal (but now adjusted) cannula-carrier plane. This 
scale, located between the PED scale 160 and the vertical 
slide 120 of the cannula-carrier, permits a rotation of the 
cannula-carrier plane around its adjusted X axis ("ML" axis) 
by a measured angular distance (ROT) . (Note that once 

25 adjusted, the originally AP, ML, and DV axes will be 

designated in quotation marks in this narrative.) The physical 
axis of rotation for this adjustment is offset from the shaft 
of the guide-tube cannula itself in the Kopf stereotactic, but 
like the PED offset this offset can be compensated by a 

30 translational adjustment later on. In effect, adjustment at 
this joint rotates the cannula-carrier plane from its initial 
orientation, with its DV axis parallel to the sagittal plane 
of the cranium, to a new orientation, with its "DV" axis 
parallel to the guide-tube cannula. This new position is 

35 specified by the direction cosines of the guide-tube cannula, 
which imply a new orientation within the cannula-carrier plane 
at an angle ROT from th Y axis (DV axis) . 
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The angle ROT can be calculat d from the known 
points at the top and the tip of the guide-tube cannula by 
considering this cannula-carrier plane (Fig. 6B) . The cannula 
is a transverse line in the cannula-carrier plane that 
constitutes the hypotenuse of a triangle of which the height 
(parallel to the DV axis) is: 



10 



15 



20 



(7) 

A - DVtop - DVtip 

and the base (the transverse line across the axial plane 
(MVAP plane or Z plane)), as calculated in (3) above is: 

C - SQRT{ APtop - APtip) 2 + (MLtop - MLtip)* ) 

The angle of the hypotenuse from the DV axis is therefore: 

ROT - (-APorient) * arctan { (SQRT(APtop-APtip) ' ♦ (MLtop - 
MLtip) 2 ) ) / (DVtop-DVtip) 



In Equation (9) , the tangent (of which ROT is the 
arctan) can obviously lie between negative and positive 

25 infinity. Por each tangent value the default arctan value is 
between -ir/2 and w/2, although there is also a distinctive 
angle beyond this range. Hote that if DVtop > DVtip, the 
tangent will always be positive, and the angle prescribed by 
Equation (9) will lie between 0 and w/2 (or between -*/2 and 

30 but the default choice is the angle between 0 and w/2) . 

As long as the guide-tube cannula 20 approaches a tangent 
point between the most posterior point and the most ventral 
point on the circle, APorient will be negative and the angle 
positive. Since the adjustment on the Kopf stereotactic moves 

35 the movable line anterior to the stationary line to set this 
angle, this sign is appropriat . If th approach is from the 
ant rior directi n, th guide-tube cannula will be tilt d in 
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the opposite mann r, APorient will be positive, and the angle 
will appropriately be negative. 

An approach in which DVtop < DVtip is not permitted 
by the Kopf stereotactic equipment, and cannot be 
5 appropriately accommodated by Equation (9) . 

4. A fourth step is to specify the "J-angle 11 , the 
angle of rotation at which the arc of curved nitinol tubing 
will emerge. 

The physical rotation of axes by the angles PED and 

10 ROT outlined in Part V Step 3 above brings the physical guide- 
tube cannula parallel to the position specified by its 
direction cosines r as provided by the Core Algorithms. The 
appropriate translational adjustments in three dimensions 
outlined in Part VI below will bring it into position as 

is specified by the Core Algorithms. Although the cannula is now 
tangent to the circle planned by the Core Algorithms, it is 
also tangent to any other circle with a point of tangency at 
this point; an entire torus of tangent circles could be formed 
around the point of tangency, and the cannula would be tangent 

20 to any of them. Thus, as the curved nitinol tube emerges from 
the fully positioned guide-tube cannula, its plane of 
curvature can be oriented in any direction; but only one 
direction is correct. Therefore the algorithms must now 
specify that the curved nitinol tubing must be set at an angle 

25 of rotation (J) from the cannula-carrier plane. This is a 
convenient plane of reference for the J-angle, because when 
the PED and ROT angles are returned to 0, it is parallel to 
the sagittal plane. 

The J-angle is calculated by considering that the 

30 plane of the circle must intersect the cannula-carrier plane, 
and that the guide-tube cannula must be included in both the 
plane of the circle and the cannula-carrier plane. As the 
curved portion is rotated (conceptually) from the cannula- 
carrier plane to the plane of the circle, its plane therefore 

35 moves through precisely the J-angle between those two planes. 
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In general, the angle 6 between any two planes 

(1) 

A 1 *X+B 1 *Y+C 1 *2+D 1 -0 

5 

and 
(2) 

A 2 *X + B 2 *Y + C 2 *Z + D 2 -0 

io is given by the formula 
(3) 

COS* « | A x * A 2 + B x * B 2 + C x * C 2 | / I (SQRTCAj 2 ♦ B x 2 + C x 2 > ) * 
(SQRT<A 2 2 ♦ B 2 2 ♦ C 2 2 )) ] 

(Pr otter et al., supra, p. 549) 

15 

For the circle plane: 
(4) 

A x - AA - H / SQRT(H* ♦ I 2 + J 2 ) 
B 1 « BB » I / SQRTtH 2 ♦ I 2 + J 2 ) 
20 C x - CC - J / SQRT (H 2 ♦ I 2 + J 2 > 

For the cannula-carrier plane, after the PED 
adjustment (with or without the ROT adjustment) , all values of 
X and Z must satisfy the relationship: 

25 

(5) 

(AlphaDirecCos / GammaDirecCos) - <X 2 - X) / <Z 2 - Z) 

where X 2 and Z 2 are the ML and AP values of a point known to 
30 be on the cannula -carrier plane. Substituting the guide-tube 
cannula tangent point for X 2 and Z 2 in Equation (5) , and 
rearranging to conform to the format of (2) , the formula of 
the cannula-carrier plane is: 

35 (6) 

(1) * X + (0) * Y - (AlphaDirecCos / GammaDirecCos) * Z - (M^g tang - 
AP gtang * (AlphaDirecCos / GammaDirecCos)) - 0 



and 
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(7) 



A 2 - 1 
B 2 « 0 



C 2 - (- AlphaDirecCos / GammaDirecCos) 

Substituting (4) and (7) in (3) provides: 



(8) 



cos (J-angle) - | AA * (1) + BB * (0) + CC * (-AlphaDirecCos / 
10 GanunaDirecCos) | / t <SQRT(AA 2 + BB 2 + CC 2 ) ) * (SQRT(U) 2 + (0) 2 ♦ 
(-AlphaDirecCos / GammaDirecCos) 2 ) ) ] 

or , simpl i f y ing : 
(9) 

15 cos (J-angle) - | AA - CC * (AlphaDirecCos / GanunaDirecCos) | / 

[ (SQRT(AA 2 + BB 2 + CC 2 ) ) * (SQRTd + (AlphaDirecCos / GammaDirecCos) 2 )) ] 



To account for a full 360° range of possible 
J-angles, a sign must be provided for this expression. 

20 Counting in order from the Anterior axis through quadrants I 
(+Z, +X: Anterior Left) and II (-Z, +X: Posterior Left), the 
cosines go from 1 to -1, and the angles provided by the cosine 
function go from 0 to 180°. Continuing through quadrants III 
(-Z, -X: Posterior Right) and IV (+Z, -X: Anterior Right), the 

25 cosines return from -1 to 1, and the angles decrease from 180° 
back down to 0°. This second series is to be given a negative 
sign, indicating that the J-angle adjustment which rotates the 
guide-tube cannula moves the center of the circle to the right 
of the cannula-carrier plane. 

30 If the center of the circle is to the right of the 

plane, then its ML will be less than the ML of a point on the 
plane at the same AP. To find such a point, the AP of the 
center of the circle (zz) is substituted for Z in the formula 
for the cannula-carrier plane (Equation (6) above) , and the 

35 equation is rearranged and simplified to give: 
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(10) 

X - MLgtang + (AlphaDirecCos / GammaDirecCos) * izz - AP gtan g) 

If xx < X, then the arc should turn to the 
subject's right and it should have a negative sign; otherwise 
it should be positive. Therefore the J-angle is given by: 

(ID 

J - SIGN{ xx - (MLgtang ♦ (AlphaDirecCos / GammaDirecCos) * 
izz - AP^^) ) ) * arccos { | AA - CC * (AlphaDirecCos / GammaDirecCos) | 
/ [ (SQRT(AA 2 + BB 2 ♦ CC 2 ) ) * (SQKT(1 + {(AlphaDirecCos / 
GammaDirecCos) ) 2 ) ) ] } 



VI. Machine-dependent output Algorithms: From the planned 
position of the guide-tube cannula and the planned 
stereotactic angular settings to the translational adjustments 
necessary to bring the guide-tube cannula into the final 
position required by the Core Algorithms. 

There are a number of ways to derive the necessary 
adjustments, depending in part on the type of stereotactic 
apparatus to be employed in surgery. The example below 
describes one technique adapted to the Kopf stereotactic 
apparatus which can be used for monkey or small-animal 
surgery. All angles have been left in radians, although in 
the software itself they must be converted into degrees to 
conform with the units used in typical stereotactic scales. 

The Kopf stereotactic apparatus provides for three 
calibrated translational adjustments. One of these (AP) is at 
the base of the cannula-carrier, where it is clamped onto the 
stereotactic arm 110. The second (DV) is on the vertical 
slide 120, and the third (ML) is on the horizontal slide 195. 
Normal stereotactic practice involves "zeroing" the apparatus. 
This is necessary because the cannula can be attached to the 
cannula-carrier in a c ntinuous rang f diff r nt positions. 
It is therefore necessary in normal st re tactic practice to 
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"find" th tip of the cannula by setting the translational 
adjustments to bring it to a known position, usually earbar 
zero. 

In the following descriptions, "ML" and "DV" will 
5 be written in quotation marks when referring to the angled 
cannula, because after the PED and ROT angular adjustments 
have been made the "ML" and "DV" adjustments made on the 
originally vertical and horizontal slides are no longer truly 
ML and DV adjustments within the stereotactic coordinate 
10 system, which is still the frame of reference for the 
algorithms. 

If the angles PED and ROT are small enough, the tip 
of the angled cannula can be zeroed as usual. (But see Part 
VIII Step 1 for the alternative case.) The settings on the 

is "ML", "DV", and AP scales at earbar zero are designated MLZ, 
DVZ, and APZ, respectively. Even though the zero points are 
known, however, the proper "ML" and "DV" adjustments that need 
to be made to move from earbar zero to the prescribed 
implantation position cannot be calculated as straightforward 

20 additions or subtractions, as is the case for a non-angled 

cannula. The formulas for calculating the final settings are 
presented below. 

1. A first step is to make the insertion-axis 
adjustment. As shown in Fig. 7, which represents a projection 

25 of the guide-tube cannula 20 onto the sagittal plane (DV/AP 
plane; YZ plane), when PED » 0, a "DV" adjustment of d 
millimeters on the vertical slide 120 of the stereotactic 
instrument moves the tip of the cannula along the insertion 
axis of the cannula from a position O (x x , y lf z 2 ) 

30 representing the origin, to a position P (x 2 , y 2 # ^2) 

representing the target point. For any such two points along 
the axis of the cannula, regardless of the value of PED, the 
DV values y x and y 2 and the distance between the points d must 
satisfy the equation: 

35 (1) 

BetaDirecCos gtang • <y 2 - Yi* I <* 
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Equation (1) implies, by substitution and 
rearrangement, that the "insertion-axis adjustment 91 that will 
bring the tip of the cannula from earbar zero (0, 0, 0) to the 
final prescribed position (MLtip, DVtip, APtip) has the 
5 length: 

(2) 

d » DVtip / BetaDirecCos gtang 

io If the reading on the "DV" scale (the vertical 

slide 120) was originally DVZ, then after the adjustment it 
will be: 

(3) 

15 AX - DVZ - (DVtip / BetaDirecCoSgt^g) 

The guide-tube cannula direction cosines are 
signed, as noted above (Part III Step 2 Equation (17)), to 
reflect the signed distance from the cannula "top" to the 

20 cannula ••tip". Correspondingly , the quantity (DVtip / 

BetaDirecCos gtang ) measures the signed distance, in the same 
direction, from the DVtip level to the zero level. But the 
adjustment desired is the AX adjustment on the vertical scale, 
which is the distance from the zero level to the DVtip level. 

25 Therefore the quantity d must be subtracted as shown. 

2. A second step is to make the cross-field 
adjustment. For the adjustment along the horizontal slide 
195, it is helpful to picture the projection of the tip of the 
30 guide-tube cannula 20 onto the axial plane (MVAP plane; XZ 
plane) as it moves during the adjustments (Fig. 8). 

Point 0 represents the location of earbar zero, and 
line DE represents the distance MLtip from earbar zero. After 
the angle adjustments FED and ROT have been made and the 
35 cannula has been zeroed, its tip is at point O. The 

ins rtion-axis adjustment (just calculated in Part VI Step 1) 
along the direction of the guide-tube cannula moves th tip f 
th cannula along the axis of the cannula d millim t rs from a 
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position (x lf y lf z x ) to a position (x 2 » Y 2 » 2 2> • ("H" 5 
projections of thes points onto the XZ plane are labeled O 
and A, respectively, in Fig. 8.) This distance d must satisfy 
the equation: 

(1) 

AlphaDirecCos gtang «= (x 2 - Xj_) / d 

This implies, by substitution and rearrangement, 
that the "DV M adjustment d along the direction of the guide- 
tube cannula that brings the tip of the cannula from earbar 
zero (0, 0, 0) to a position (ML, y, z) has the length: 

(2) 

d - ML / AlpbaDirecCos gtang 

and that the effect of such an adjustment on the ML position 
of the tip will be: 

(3) 

ML - d * AlphaDirecCos gtang 

The line segments OA, OB, and OC all cover the same 
X distance (ML - d * AlphaDirecCos gtang ) in the projection 
diagram (Fig. 8) ; for line OB along the X axis itself this 
distance is equivalent to its length. If this X distance were 
the X distance prescribed by the Core Algorithms as MLtip, no 
further adjustment would be necessary. If some other X 
distance is prescribed, however, then a further adjustment CD 
will need to be made along the horizontal slide 195 (the 
"cross-field" slide) of the stereotactic instrument. The 
projection of this adjustment on the XZ plane (line segment 
CD, and by extension line segment OD) is perpendicular to the 
projection of the insertion-axis adjustment AO, because the 
horizontal slide 195 is perpendicular to the vertical slide 
120 on the st reotactic instrument (Figs. 4-5). Therefore, 
the angle PED is repeated in the triangl OBC, f rmed by th 
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inters ction between OD and the extension of line segm nt AB, 
as illustrated. Therefore the line OC has the length: 

(4) 

5 OC - (d * AlphaDirecCofl gtang ) / ( cos (PBD) ) 

This line segment, which marks the effect of the 
insertion-axis adjustment on the X position, and the line 
segment CD, which marks the additional distance to be moved 
10 along the horizontal slide to reach the final X position, 

compose the line OD in triangle ODE. The angle PED is again 
repeated in triangle ODE, so that: 

(5) 

15 COS (PBD) - DB / OD 

The known lengths of DE (MLtip) and OC (Equation 
4), together with the unknown length CD can be substituted in 
Equation (5) to provide: 

20 

(6) 

cos (PBD> - MLtip / [ CD + (d * AlphaDirecCos gtang ) / (cos (PBD) ) ] 
Simplifying and rearranging Equation (6) provides: 

25 

(7) 

CD * cos (PBD) - MLtip - (d * AlphaDirecCoSgt^jg) 
which provides the value of CD: 

30 

(8) 

CD - [ MLtip - (d * AlphaDirecCos gtang ) ] / [ cos (PBD) ] 

The value of d can now be substituted from Part VI 
35 Step 1 Equation (2) above, to yield: 
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(9) 

CD - t MLtip - { (DVtip / BetaDirecCos gtang ) * AlphaDirecCos gtang ) ] 
/ [cos (PBD) ] 

5 In practice, to meet the requirements of the sign 

convention (right negative, left positive) , and to permit use 
of the carrier on either arm of the stereotactic, this formula 
can be modified by using the variable SideArm (introduced in 
Part V Step 3 Equation (5)): 

10 

(10) 

CD - SideArm * [ MLtip - < (DVtip / BetaDirecCos gtang ) * 
AlphaDirecCoB gtang ) ] / [ cos (PBD) ] 

15 if the reading on the "ML" scale (the horizontal 

slide) was originally MLZ, then after the adjustment the final 
reading will be: 

(11) 

20 CR - MLZ + CD 

or: 
(12) 

CR - MLZ + SideArm * [ MLtip - ( (DVtip / BetaDirecCos gtang ) * 
25 AlphaDirecCos gtang ) ) / t cob (PBD) ] 

3. A third step is to make the along-field 
adjustment. This involves positioning the carrier 
(essentially consisting of items 120, 160, 170, 190, 195) 
30 along the stereotactic arm 110. For the final adjustment, it 
is necessary to take into account the effects of the 
insertion-axis adjustment and the cross-field adjustment on 
the along-field position of the guide-tube cannula tip. 

The insertion-axis adjustment described above in 
35 Part VI Step 1, of magnitude d along the direction of the 
guide-tube cannula, moves the tip of the cannula along the 
axis of the cannula d millimeters from a position (x 1# y x , z x ) 
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to a position (x 2 , y 2 , * 2 ) • Tnis distance d must satisfy th 
equation: 

(1) 

GammaDirecC08g Can g - (z 2 - *i> / d 

This implies, by substitution and rearrangement, 
that the insertion-axis adjustment d along the direction of 
the guide-tube cannula, which brings the tip of the cannula 
from earbar zero (0, 0, 0) to a position (X, Y, AP) , has the 
length (dashed line PQ in Fig. 6): 

(2) 

d - AP / GammaDirecCoSg tan g 

and that the effect of such an adjustment on the AP position 
of the tip will be: 

(3) 

AP - d * GamnaOirecCo8g tang 

The effect of the cross-field adjustment (dashed 
line CP in Fig. 7) can be calculated directly from triangle 
CDF, since: 

(4) 

sin (PBD) - CP / CD 
so that 
(5) 

CP - CD * ( sin (PBD) ) 



or, substituting for CD from Part VI Step 2 Equation (10) 
above: 
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(6) 

CF = {SideArm * [ MLtip - ( (DVtip / BetaDirecCos gtang ) * 
AlphaDirecCos gtang ) ) / [ cos (PED> ] } * { sin (PBD) ) 

The tip of the cannula must be placed at a distance 
APtip from earbar zero as prescribed by the Core Algorithms. 
The adjustments AP and CF (Equations (3) and (6)) have already 
been made, so the remaining distance to be covered is: 
(7) 

APslide - APtip - AP - CP 

(8) 

APslide - APtip - (DVtip / BetaDirecCos gtang ) * GammaDirecCos gtang 
- { SideArm * [ MLtip - < (DVtip / BetaDirecCos gtang > * 
AlphaDirecCos gtang ) ] / [ cos (PED) ] } * ( sin (PBD) ) 

If the reading on the AP scale (on the stereotactic 
arm) was originally APZ, then after the adjustment the final 
reading will be: 
(9) 

APadj - APZ + APslide 

or: 
(10) 

APadj - APZ + APtip - (DVtip / BetaDirecCos gtang ) * 
GamroaDirecCos gtang - { SideArm * [ MLtip - ( (DVtip / BetaDirecCos gtang ) * 
AlphaDirecCos gtang ) ] / [ cos (PBD) 3 } * (sin (PED) ) 



VII. Machine-dependent output Algorithms: From a series of 
trajectory points and characteristics of the implantation tool 
to prescribed tubing measurements. 

1. A first step is to determine, for the curved 
part of the trajectory, the length of the arc extending from 
one of the parallel planes to the next. Again, there are a 
number of ways to do this, as there are for doing most of the 
tasks to follow. One useful way is illustrated below. 
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If the points are d signated (ML^, DV i-lf AP^) and (ML ir 
DV if APi) , th length of the chord b tween the two points 
(Fig. 9) is: 



5 (1) 

CHORD £ - SQRT ( (ML^ - ML^) 2 + i^i ' ^^l 2 + (AP ± - AP^) 2 ) ) 

The line bisecting the angle between the radii to 
the two points also bisects the chord and forms a right angle 
10 with it. The hypotenuse of the right triangle on either side 
of the bisector is the radius (rad) of the circle, and the 
base is half the chord. Therefore the sine of the angle 
(theta / 2) between radius and bisector can be written as the 
base divided by the radius. 

15 

8in (theta / 2) - (SQRT ( (HL ± - ML^) 2 + (DVi - DV^) 2 + ihP L - 
AP^) 2 ) ) / 2) / rad 
or 



20 (2) 

sin (theta / 2) - (SQRT ( (ML A - ML^) 2 + (D^ - DV^) 2 ♦ (APi - 
AP^) 2 )) ) / (2 * rad) 

Twice this angle is the angle that subtends the arc 
25 and therefore the angular length of the arc itself in radians. 



(3) 

theta - 2 * arcsin((SQRT ( (ML^ 
AP i-l )2)) > / < 2 * rad,) 

30 

This angular length 
linear length of the arc. 



- ML^) 2 + (DV ± - D^-i) 2 + (AP ± - 

times the radius gives the 



(4) 

35 ARCi - 2 * rad * arcsin((SQRT ( (ML^ - ML^) 2 + iW ± - DV^) 2 + (AP A 

- AP i-l >2)) > 1 (2 * rad)) 



2. A second st p is t cumulate thes arc lengths 
all along the curv d trajectory, by adding the arc lengths 
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calculated above. The resulting series of increasing lengths 
provid s the basis for calculating the arc length from any one 
point to any other point, simply by subtracting one cumulated 
arc length from the other. Two such sets are constructed, for 
different purposes, one beginning with the most posterior 
point on the circle, and the other beginning with the guide- 
tube cannula tangent point. 

3. A third step is to determine, for the straight 
parts of the trajectory, the length of the line extending from 
the tangent point to the intersection point with each of the 
parallel planes. This provides for each tangent a series of 
lengths that is comparable to the series for the arcuate 
portion derived above. 

If the intersection point is designated (ML i# DV A , 
AP L ) , the distance along the guide-tube cannula, for example, 
is simply: 

(1) 

LBG ± - SQRTUAPi - AP gtang ) 2 + (MI^ - MI^ tang > 2 ♦ (DV ± - DV gtang > 2 ) 

4. A fourth step is the calculation of the points 
at the most extreme Z axis levels represented on the circle. 
A number of lengths are required as specifications, both for 
the metallic elements of the implantation tool and for the 
flexible implantable tubing itself. One additional 
calculation (this step) and several additional entries (next 
step) must be made to accomplish this. 

There are a number of ways to visualize and 
calculate the most posterior and anterior points on the 
circle. One efficient way is illustrated below. 

First, consider the circle plane (Part II Step 4 
Equation 1) : 

(1) 

AA*X+BB*Y+CC*Z+DD-0 
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Any plane parallel to the circle plane will have 
coefficients in the same proportion as coefficients (AA, BB, 
CC) , which represent the direction components of any line 
perpendicular to the plane (Protter et al., supra, p. 546). 
The direction cosines of any such line can be computed from 
these numbers as: 

- AA / SQRT(AA 2 + BB 2 + CC 2 ) 

- BB / SQRT (AA 2 + BB 2 + CC 2 ) 

- CC / SQRT ( AA 2 ♦ BB 2 + CC 2 ) 

In particular, with reference to Fig. 10 , consider 
the plane through the origin of the coordinate system and 
parallel to the circle plane. Translate the circle to this 
plane, with its center at the origin (O) . The posterior point 
on the circle (C) , at AP level P, will be in the same plane 
that contains the Z axis and the normal axis of the circle 
(the line normal to its plane and passing through its center) . 
This plane is pictured in Pig. 10. Since the normal to the 
circle is in the plane, the circle necessarily appears as a 
straight line (CR, its projection onto this plane) . The line 
segment OC is a radius (rad) of the circle. In the triangle 
OCP, it is evident that: 

(5) 

cos (»/2 - Gamma) - OP/rad 

so that the AP coordinate of the furthest posterior point (C) 
on the circle is: 

(6) 

OP - rad * cos ( ir/2 - arccos (GaiwnaDirecCoBjjQpM^) ) 
Substituting from Equation (4) above pr vides: 



(2) 

AlphaDirecCoSfloaMftL 

(3) 

BetaDirecCos N0RMAL 

(4) 

GairanaDirecCoSflQjgjkL 
(Love, supra, p. 128) 
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(7) 

OP - rad * cos ( w/2 - arccos (CC / SQRT(AA 2 + BB 2 + CC 2 )) ) 

To calculate the most anterior and posterior planes 
5 containing a point on the circle, this length can be 

(respectively) added to or subtracted from the Z value of the 
actual center of the circle: 

(8) 

10 AnteriorAP - zz + rad * cos ( w/2 - arccos (CC/SQRT{AA 2 + BB 2 + CC 2 )) ) 

Poster iorAP - zz - rad * cos { ir/2 - arccos (CC/SQRT(AA 2 + BB 2 + CC 2 ) ) ) 

15 (10) The ML and DV values for these points (AnteriorML and 
AnteriorDV; PosteriorML and PosteriorDV) are available by 
application of Part IX Step 4 Equation (10) and Part II Step 6 
Equation (2). 

20 5. A fifth step is the entry of certain additional 

values by the user. 

a. Selection of proximal and distal points defined 
as the limits of the "target structure* 1 along the trajectory. 

In Part III Step 1 the points of tangency of the 

25 guide-tube cannula and the wire were selected by the user. 
These points condition the course of the trajectory itself, 
and were therefore included in the Core Algorithms above. The 
points of tangency do not necessarily — and in practice usually 
will not — represent the limits of the structure to be 

30 accessed. Once the trajectory is established, an independent 
pair of points on that trajectory has to be selected by the 
user to represent the beginning and end of the target 
structure. 

As above, there are a number of ways in which the 
35 two points might be designated. One useful way is for the 

user simply to s lect n of the parallel planes for which the 
traject ry^ int rsection points have been calculated. One f 
these points is designated as (InitStructML, InitStructDV, 
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InitStructAP) , and the other as (FinalStructML, FinalStructDV, 
FinalStructAP) indicating their X, Y, and Z coordinates, 

respectively. 

b. Selection of the entry and exit points on the 

body. 

Another pair of points to be designated, once the 
trajectory is established, are the body surface point at which 
the guide-tube cannula enters the body, and the exit point at 
which the wire emerges from the body. One way to provide 
these points would be to perform an image analysis of the 
scanner images, to permit the analysis to check after 
calculating each intersection of the trajectory with a coronal 
plane, to determine whether that point lies beyond the 
cranium. A simpler way is for the user to select the proper 
planes from the GUI output described above. One of these 
points is designated as (EntryHL, EntryDV, EntryAP) , and the 
other as (ExitML, ExitDV, ExitAP) indicating their X, Y, and Z 
coordinates, respectively. 

For convenience, the AP levels of the labeled 
points defined on the circle and along the trajectory are 
listed here: 

Parallel plane nearest the guide-tube cannula entrance to the 
cranium: 

~ EntryAP 

AP coordinate of the most posterior point on the circle: 
« PosteriorAP 

AP coordinate of the tangent point at which the guide-tube 
cannula contacts the circle: 
«= APgtang 

Parallel plane nearest the trajectory entrance to the target 
structure: 

= InitStructAP 
Parallel plane nearest the trajectory exit from the target 
structure: 

- FinalStructAP 
AP co rdinate of the tangent point at which the wire leaves 

the circle: 
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= APwtang 

AP coordinate of the most anterior point on the circle: 
- AnteriorAP 

Parallel plane nearest the guide-tube cannula exit from the 
5 cranium: 

= ExitAP 

c. The following values as defined below are also 
to be entered by the user. 

The desired distance, when the guide-tube cannula and curved 
10 nitinol tubing are fully advanced, between the lowest collar 

mount (bearing the cannula) and the second collar mount 

(bearing the curved nitinol tubing): 
= FirstToSecondDisk 

The desired length of wire to be showing above the cranium 
15 after full advancement: 

» WireAboveHead 

Normally the following value should be set to zero, but if 

APwtang •goals AnteriorAP or Poster iorAP, then the point at 

which the wire emerges from the cranium cannot be calculated 
20 automatically and must be estimated by the user: The 

estimated length of nitinol wire from its point of tangency to 

the cranium: 

-WirelnHead 

Similarly, if APwtang nearly equals AnteriorAP or PosteriorAP 
25 (nearer than SliceThickness) , then it is convenient to include 
among the sample coronal planes an extra plane containing the 
point at which the wire reaches the surface of the cranium. 
For this purpose, the user may provide the distance from such 
an extra plane to AnteriorAP or PosteriorAP: 
30 - SpecialWirelncrement 

Likewise, if APgtang nearly equals AnteriorAP or PosteriorAP 
(nearer than SliceThickness) , then it is convenient to include 
among the sample coronal planes an extra plane containing the 
point at which the cannula reaches the surface of the cranium. 
35 For this purpose, the user may provide the distance from such 
an extra plane to AnteriorAP or Poster i rAP: 
= SpecialGuid Increment 
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6. A sixth step is the summation of values to 
provide lengths for th preparation of tubing to meet the 
requirements of the Core Algorithms. 

The designated coronal planes (Part VTI step 5b) 
are used, in conjunction with the arc lengths and leg lengths 
determined in Part VTI Steps 1-3, to summate the arc length 
and/or trajectory length from one designated point to another. 
For making these calculations, reference is made to two sets 
of cumulative lengths. One set (ARC) is defined, for planes 
from one point of tangency to the other, as the cumulative 
length of the arc of the circular path, starting from the 
point of tangency with the guide-tube cannula and ending with 
the designated point; for planes beyond the points of 
tangency, this set is defined as zero. The second set (LEG) 
is defined as zero for all points along the trajectory between 
the two points of tangency, inclusive? for points along the 
trajectories of the guide-tube cannula or the emergent wire, 
this set is defined as the distance from the given point to 
the respective point of tangency. 
a. Circle: 

(As written, the following definitions assume the 
target is entirely included in the circular portion of the 
trajectory; the first three lengths are therefore useful only 
if this is the case. In each case, what is meant by the 
definition is the difference between the cumulative arc length 
along the circle at the later point and the cumulative arc 
length at the earlier point.) 

Posterior to targeted structure: 

■» ARC InitStruct - ARCposterior 

Within target: 

= ARC pinalstruct - ARC InitStruct 

Target to anterior point: 

- ARCantertor - ARC pinalStrucC 
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Posterior to Anteri r: 

b. Trajectory: 

(The following definitions are more general, making 
no assumption about the overlap of the target structure and 
the circular portion of the trajectory. If the trajectory is 
still within the target beyond the point of wire tangency, for 
example, the "TargetToWireTangency" becomes a negative number 
showing how far the target reaches beyond the point of wire 
tangency.) 

Cran iumToTangency : 

« LEG Bntry - LEG gtang 

TangencyToTarget : 

m ARC InitStruct - ARC gtang - ^GmitStruct 

WithinTarget: 

~ ARC pinalStruct - ARCmitStruct + 
LE6 PinalStruct + LEG InitStruct 

TargetToWireTangency : 

= ^RCwtang " ARC PinalStruct ~ ^^FinalStruct 

WireToCranium : 

« LEGBxit " ^wtang 

TangentToTangent : 

" ARC wtang " ^Ggtang 

c. Compound implantable tubing: 

Combinations of the Trajectory lengths (Part VII 
step 6b) are used in the construction of the compound 
implantabl tubing, as follows. 
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The n l ading" 1 ngth of non-pat nt PAN-PVC tubing to be pulled 
in and implanted (on the guide-tube cannula "sid " of the 
target structure) by retraction of the implantation tool: 
Posterior Occlude - CraniumToTangency + TangencyToTarget 

5 

The length of patent PAN-PVC to be pulled in and implanted (in 
the target structure) by retraction of the implantation tool: 
Patent Span - withinTarget 

10 The trailing length of non-patent PAN-PVC tubing to be pulled 
in and implanted (on the emergent wire "side" of the target 
structure) by retraction of the implantation tool: 

Anterior Occlude - TargetToWireTangency + WireToCranium 

15 d* Implantation tool: 

Specifications for the implantation tool 10 itself 
are required for manufacturing. 



20 



25 



Guide-tube cannula: 

Length of guide-tube cannula 20 in the block 15 and the lowest 
collar mount: 

GuideCannulalnBlockAndLowestDisk « Mount ingBlockThicknese + 
CollarMountThickness 



Length of guide-tube cannula 20 shoving, from the bottom of 
the block 15 to the top of the cranium, when the tip of guide- 
tube cannula is in final position: 

GuidePromBlockToCranium - LB - SQRT( (APtip - BntryAP) 2 + (MLtip - 
30 KntryML) 2 + (DVtip - BntryDV) 2 ) 

Length of guide-tube cannula 20 inside the cranium 150, from 
cranium to tip, in final position: 

Guide Ins ideCr an ium . SQRT( (APtip - BntryAP) 2 + (MLtip - BntryML) 2 + 
35 (DVtip - BntryDV) 2 ) 



Total guide-tube cannula 20: 

* Guid CannulalnfilockAndLowestDisk + GuidePromBlockToCranium + 
GuidelnsideCranium 



WO 96/03087 



PCT/US95/09391 



82 

Nitinol tubing: 

Tubing in the block 15 and lower collar mount 16 and middle 

collar mount 14: 

NitinolTubinglnBlockAndDisks - Mount ingBlockThickness + (2 * 
CollarMountThickness) 

Tubing from the bottom of the block 15 to the tip of the 
guide-tube cannula 20: 

TubinglnGuide - LB - CurvelnGuideLength + CurvelnGuideArcLength 

Length of the nitinol tubing 28 from the tip of guide-tube 
cannula 20 to the emergence of nitinol wire 24 from the 

nitinol tubing 28: 

TubingTravel - TangentToTangent - CurvelnGuideArcLength + 
StraightWireArcLength 

Total nitinol tubing 28: 

- NitinolTubinglnBlockAndDisks ♦ TubinglnGuide + TubingTravel 

The nitinol tubing 28 should be manufactured with the 
curvature extended slightly beyond what is necessary — the 
slight excess includes the concealed CurvelnGuideArcI^ngth, 
and a little extra which will be concealed in the guide-tube 
cannula 20 and have no influence on the trajectory: 
PlannedBxcessArc - 1 + INT (CurvelnGuideArcLength + 0.5) 

Nitinol tubing straight part: 

- NitinolTubinglnBlockAndDisks + TubinglnGuide - PlannedBxcessArc 

Nitinol tubing curved part: 

» TubingTravel + PlannedBxcessArc 

Nitinol wire 24: 

The length of nitinol wire 24 in the block 15 and the collar 
mounts 16, 14, 12, plus the length of the stop 30 between the 
lower and middl c liar m unts 16, 14 when the wire is fully 
advanced: 
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NitinolWirelnDisks » Mount ingBlockThickness + (3 * 
CollarMountThickness) + FirstToSecondDisk 

The length of nitinol wire 24 from the bottom of the block 15 
5 to the tip of the nitinol tubing 28: 

WirelnTubing - TubinglnGuide + TubingTravel - StraightWireArcLength 

The length of nitinol wire 24, from the nitinol tubing 28 to 
the point of emergence from cranium, plus the amount to be 
10 shoving above cranium: 

WireTravel - WireAboveHead + WireToCranium 

or, if WirelnHead is given a non-zero value: 
WireTravel - WireAboveHead + WirelnHead 

15 Total nitinol wire 24: 

- NitinolWirelnDisks + WirelnTubing + WireTravel 

Guide-rail: 

20 Third guide-rail (gauge 22): 

- PirstToSecondOisk + 2 * CollarMountThickness + WireTravel 



VIII. Alternative input methods 

25 

1. A first alternative is an alternative method 
for zeroing the guide-tube cannula 20. 

Mote that in Part VI it was assumed that the 
stereotactic carrier could be zeroed in the normal way, at 

30 earbar zero. For placements that involve relatively large 
angles PED and ROT, however, it is not possible or not 
convenient to reach both earbar zero and the point (MLtip, 
DVtip, APtip) at which the guide-tube cannula is to be set for 
the implantation. In this case, it is necessary to set up a 

35 second, non-angled "marker" carrier which can be zeroed at 

earbar zero and then conveyed to (MLtip, DVtip, APtip) . The 
guide-tube cannula is then positioned t t uch the same point, 
and the settings on the horiz ntal slide (CRz ro) , vertical 
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slide (AXzero) , and stereotactic arm (APzero) are recorded 
directly. 

It is helpful, for purposes of checking the 
settings and the formulas, to back-calculate the location of 
5 earbar zero in this case, but it is not necessary. This back- 
calculation requires a number of formulas, the derivation of 
which can be shown by considerations similar to those in Part 
VI above. Therefore the formulas leading to this back- 
calculation are merely presented below, without derivation, 
xo First, it is useful to express the effects, on X f 

Y, and Z position (all else equal), of an insertion-axis 
adjustment that increases "DV" position by 1 millimeter: 

(1) 

15 AXmlEff ect - - AlphaDirecCo8 gt:an g 

AXdvEffect - - BetaDirecCofig^ng 
AXapBf feet - - GammaDirecCos gtan g 

Likewise, it is useful to express the effects, on X, Y, and Z 
20 position (all else equal) , of a cross-field adjustment that 
increases "ML" position by 1 millimeter: 



(2) 

CRmlBffect - SideArm * cos(PBD) 
25 CRdvBffect - 0 

CRapEffect - sin(PED) 



Since the stereotactic arm remains parallel to the 
Z axis, an along-field adjustment of 1 mm (all else equal) 
does not change ML or DV position, and changes AP position 
simply by the nominal 1 mm. 

If adjustments are to be made on all three scales, 
then "all else" is not "equal". Since there is only one way 
to adjust DV, this adjustment is made first, conceptually. 
The insertion-axis adjustment that will bring the guide-tube 
cannula tip from DVtip to the DV of earbar zero is: 
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(3) 

AXtoZero - - DVtip / AXdvBf feet 

The effects of this move in the other two 
5 directions are: 

(4) 

AXtoZeroMLef feet - AXtoZero * AXmlBf feet 

10 (5) 

AXtoZeroAPeffect - AXtoZero * AXapBf feet 

Since there is only one remaining way to adjust ML, this 
adjustment is made next, conceptually. The cross-field 
is adjustment that will bring the guide- tube cannula tip from its 
new position to the ML of earbar zero is: 

(6) 

CRtoZero - {- NLtip - AXtoZeroMLef feet) / CRmlBffect 

20 

The effect of this move in the AP direction is: 

(7) 

CRtoZeroAPeffect - CRtoZero * CRapEffect 

25 

The along-field adjustment that will bring the 
guide-tube cannula tip from its new position to the AP of 
earbar zero is: 

30 (8) 

APadjtoZero - - APtip - AXtoZeroAPeffect - CRtoZeroAPeffect 

The calculated values for earbar zero are then: 

(9) 

35 MLZ - CRzero + CRtoZero 

DVZ - AXatero + AXtoZero 
APZ - APaero + AP djtoZero 
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If the values of Equation (9) are used as MLZ, DVZ, 
and APZ to calculate the AX, CR, and APadj values , th y will 
arrive at the original inputs made at the beginning of this 
step. 

5 

2. A second alternative is an alternative input 
method, designed for making inputs from measurements in a 
brain atlas. The atlas may be a conventional atlas, or it may 
be a set of scanner images for which the .TIF files are not 
10 available. 

This input method allows readings to be entered 
directly from the scanner images or atlas and provides scaling 
and offsetting to convert them to the surgical stereotactic 
coordinate system with origin at earbar zero. This method 

15 (assuming the pictures available are coronal views) requires 
from the user (1) a pair of distances, one number 
(ScaleMeasured) representing a length measured on the scanner 
image or atlas, and another (ScaleReal) representing the Known 
corresponding distance in real millimeters; (2) for each of 

20 the three target points, its measured distance from the 

midline (HLmeas) and from the bottom of the picture (DVmeas) , 
and the number of its coronal section (APmeas) , in the series 
of coronal sections (in typical MRI sections this represents 
real millimeters from machine zero, but machine zero does not 

25 fall at the earbar zero AP) ; (3) the measured distance from 
the bottom of the picture for earbar zero (DVEZ) on the 
scanner image or atlas, and the number of the coronal section 
containing earbar zero (APE2) , in the series of coronal 
sections. The algorithms below convert these values into 

30 terms of real millimeters from earbar zero for the final 
entries (MI*, DV A , AP A ) , (MI*, DV B , AP B ) , (ML^, DV C , AP C ) 
suitable for the Core Algorithms. 

The pair of distances provides a scale factor for 
converting the measured ML and DV values to real millimeters: 

35 (1) 

MRItoReal - ScaleReal / ScaleMeasured 
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Th seal factor and the measured coordinates and 
earbar zero valu s are used to make the conversions for all 
three points as shown here for target point 1, and the results 
are provided to the Core Algorithms as the three original 
5 target points: 

(2) 

ML^ - [ MLmeas A ] * MRItoREAL 

(3) 

10 DV A - [ DVmea8 A ] * MRItoREAL - DVBZ * MRItoREAL 

(4) 

AP A « [ APmeas A ] - APEZ 



Analogous conversions can be made in the opposite 
15 direction to provide the coordinates of the trajectory points 
on the scanner image or atlas, for each coronal section. 



3. A third alternative is similar to the second, 
but instead of three points it takes input in the form of one 

20 point, the radius, and one angle of orientation (supplying the 
second angle of rotation by default) . 

This input method provides for a circle plane 
orthogonal to the axial plane (XZ: AP/ML) , which can be 
" targeted" on a single point (target point 2) , but placed 

25 transversely at a specified angle, rotated on its orthogonal 
(DV) axis from a conceptual "starting position" parallel to 
coronal plane. This method (also assuming the pictures 
available are coronal views) requires from the user (1) a pair 
of distances, one number (ScaleMeasured) representing a length 

30 measured on the scanner image or atlas, and another 

(ScaleReal) representing the known corresponding distance in 
real millimeters; (2) for only one target point, its measured 
distance from the midline (MLmeas B ) and from the bottom of the 
picture (DVmeas B ) , and the number of its coronal section 

35 (APmeas B ) in the series of coronal sections; (3) the measured 
distance from the b ttom of the pictur for earbar zer (DVEZ) 
on th scanner imag or atlas, and the number of the c ronal 
section containing earbar zero (APEZ) in th series of cor nal 
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sections; (4) the desired angle of rotation from the cor nal 
orientation (AxialRotation) ; and (5) the desired radius of the 
circle (WireRadius) . The algorithms below convert these 
values into terms of real millimeters from earbar zero for the 
final entries (ML A , DV A , AP A ) , (MI*, DV B , AP B ) , (MLq, DV c , AP c ) 
suitable for the Core Algorithms. 

The pair of distances provide a scale factor for 
converting the measured ML and DV values to real millimeters, 
as above, but also a scale factor for converting real 
millimeters to the scale represented by the measured values: 

(1) 

RealtoMRI - ScaleMeasured / ScaleReal 

The scale factors, the measured coordinates and 
earbar zero values, and the selected angle and radius are then 
used to make the conversions for target points 1 and 3, while 
target point 2 is simply converted as shown under Part VIII 
Step 2 above: 

(2) 

ML A - [ MLmeas B - SIGN (MUneae B ) * WireRadius * RBALtoMRI * 
cos (AxialRotation) ] * MRItoREAL 

(3) 

DV A - [ DVmeas B ♦ WireRadius * REALtoMRI ] * MRItoREAL - DVBZ * 
MRItORBAL 

(4) 

AP A - [ APmeas B ♦ WireRadius * sin (AxialRotation) 1 - APBZ 

(5) 

MI^ « [ MLmeas B ) * MRItoREAL 

(6) 

DV B « [ DVmeas B ] * MRItoREAL - DVBZ * MRItoREAL 

(7) 

AP B - [ APmeas B ] - APEZ 
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(8) 

MLc - [ MLmeas B + SIGN (MLmeaSgj * WireRadius * REALtoMRI * 
cos(AxialRotation) ] * MRItoRBAL 

5 (9) 

DV C - [ DVmeas B + WireRadius * REALtoMRI 1 * MRItoRBAL - DVBZ * 
MRItoRBAL 

(10) 

10 AP C - [ APmeas B - WireRadius * sin (AxialRotation) ] - APBZ 

The results are provided to the Core Algorithms as 
original target points 1, 2, and 3. 

15 4. A fourth alternative is an alternative input 

method, similar to Part VIII Step 3 above, designed for 
targeting a single given point, whether it is input from the 
GUI interface, from scanner images, from an atlas, or as a 
real value. It enables rotating the circle through two angles 

20 and specifying its radius. This input method is useful for 
targeting a small site with structures nearby that must be 
avoided by the entire trajectory. 

This input method assumes that the input point has 
been pre-scaled by other algorithms to expression in terms of 

25 the stereotactic coordinate system with its origin at earbar 

zero; it requires from the user (1) inputs to provide a single 
target point (MI*, DV B , AP B ) ; (2) an angle HOR through which 
to twist the circle conceptually, with that target point as 
the base, from a coronal starting position; (3) an angle SAG 

30 to lever the circle down, with the "target point" as the only 
fixed point and with the circle's angular relationship to the 
sagittal plane maintained; (4) the desired size of the radius 
(siz) ; (5) the side to be treated (SideTreated: left +1, right 
-1) . The algorithms below use these values to select the 

35 other two points (ML*, DV A , AP A ) and (MI^, DV C , AP C ) that will 
force the Core Algorithms to provide a circle that has the 
provided features. 
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The conceptual starting point is (1) a vertical 
guide-tube cannula; (2) a circle of the desired radius (siz), 
in a coronal plane, normal to the AP axis; (3) the target 
point (MLq, DV b , AP b ) fixed at the bottom center of the 
5 circle; (4) the other two points (ML A , DV A , AP A ) and (ML^, DV C , 
AP C ) originally 90° from it on either side of the circle. 
Before rotating, then, 



(1) 










ML A - 




Side Treated * siz 


(2) 










ML C - 


MLe + 


Side Treated * siz 


(3) 










DV A - 


DV C 0 


DV B + siz 


(4) 










AP A - 


AP C - 





The formulas perform the following conceptual 

steps: 

20 First, rotate the circle by HOR around its vertical 

diameter, which is parallel to DVaxis, such that the medial 
90° point moves posterior. For AP, this subtracts (siz * 
sin (HOR) ) from the medial side, and adds (siz * sin (HOR) ) to 
the lateral side. For ML, this adds ((Side Treated) * (siz - 

25 siz * cos (HOR)) to the medial point and subtracts (SideTreated 
* (siz - siz * cos (HOR)) from the lateral point. The axis 
remains vertical. The following formulas carry out this first 
maneuver: 

30 (5) 

ML A • MLjj - Side Treated * siz * cos (HOR) 

(6) 

ML C - MLg + Side Treated * siz * cos (HOR) 

35 Second, lever this vertical axis down within the 

sagittal plane by the angle SAG. The previously vertical axis 
is m ving lik the spoke of a wheel with the c nt r at th 
bottom p int. The circle moves along with it. For levering 
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the top point of the axis back and down, a point just lateral 
to the c nter com s up a little, and a point just medial comes 
down a little to lie lover than the original DV of the bottom 
point. 

5 The following intermediates simplify the formulas 

that follow: 

(7) 

angmed ■ SAG + arctan (sin (HOR) ) 

10 (8) 

anglat • SAG - arctan (sin (HOR) ) 

(») 

mul - SQRT(siz 2 + (siz * sin (HOR)) 2 ) 
15 The following formulas carry out the second maneuver: 



(10) 










- DV B 


+ mul 


* cos (angmed) 


(11) 










- AP B 


- mul 


* sin (angmed) 


(12) 








DV C 


- DV B 


+ mul 


* cos (anglat) 


(13) 










- ap b 


- mul 


* sin (anglat) 



25 

The starting point 2 (MI*, DV B , AP B ) , which has not 
moved, and the adjusted target points 1 (MI*, DV A , AP A ) and 3 
(MLc, DV C , AP C ) from Equations (5), (6), and (10) through (13) 
above, are then used as the three original target points for 
30 the Core Algorithms* 

IX. Fitting a Three-Dimensional Helical Curve to a Set of 
Points: From an ordered series of four or more points to a 
35 helical curve fit to pass through or near these points. 

Similar methods may be used to determine helical 
paths that includ curvature in thre dim nsions. The m thods 
outlined below describe th fitting of a helix to a set f 



W 96/03087 



PCT/US9S/0939r 



92 

input points, analogous to the determination of a circle from 
three points (Part II, Steps 1-3). 

1. Strategy: An Overview* 

(0) Consider n points in space, x 2 = (x x , y x , 
5 Sx),...,^, - (Xn, y n , z n ), called "the data." 

(1) Fit any curve G to or through the data, for 
example, segments of straight lines joining consecutive 
points, or a cubic spline through the data. 

(2) Parameterize the curve G with its arclength. 
10 At all data points, estimate with finite differences or 

calculate with derivatives the following objects: 

(2.1) the tangent vectors T:= dG/ds, T 1# ...,T n , so 
that they all have unit length, Bt ± | = 1; 

(2.2) the derivatives dT^ds and the curvatures 
is i^sld^/dsH ; the average curvature «: = (l/n) Ej^*! serves as 

the curvature of the helix; 

(2.3) the normal vectors N lr ...,N n , defined by 
$i : - ( V II dTi/ds I ) dT ± /ds ; 

(2.4) the binormal vectors B lf . . . ,3^, defined by 
20 the cross product B i :=T i x N ir 

(2.5) the derivatives dNj/ds and the torsions 
given by the dot product r i : = (dN i /ds) •B i ; the average torsion 
r:*(l/n)L ^_^r i serves as the torsion of the helix. 

(3) Fit a plane P to the normal vectors N lf ...,N n , 
25 for instance, by the method of total least squares (See, van 

Huff el et ai., Tfre Total Least Squares Problem: Computational 

Aspects and Analysis . Society for Industrial and Applied 
Mathematics, Philadelphia, PA, 1991). (A cheap alternative 
might use the average of the binormal vectors, (l/nJE^^ B i# 
30 as the normal to the plane, but that might prove unstable.) 

(4) Compute the projections P lr ...,P n of the data 
on the plane P. 

(5) Fit a circle to the projections in the plane 
P, for instance, by the method of Noura & Kitney f Computer 

35 Phvsics Co mmunications. 64, pp. 57-63, 1991). The radius r 0 
f the fitt d circle s rv s as th radius f th helix, and 
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the axis f the helix now passes through th center of the 
fitted circle, p rpendicularly to the plane P. 

(6) Change to coordinates (u,v,v), with (u,v) in 
the plane P and w measured along the axis of the helix. Each 

5 data point now also has coordinates (u if v i ,w i ). 

(7) Cut a piece of the helix of height w n - w x , 
and shift it parallelly to the axis, so that its endpoints lie 
at the same heights as x x and respectively. Then rotate 
the piece of helix about the axis, so that the endpoints lie 

10 at the same minimum distance from x x and x n respectively. 

(8) If desired, refine the fitted helix 
iteratively. 



2. Introduction: An Example of Subroutines With 
is Segments Between Data Points: The algorithm just outlined 

depends upon the initial curve G fitted to the data. The 
following subroutines illustrate the case with 6 consisting of 
segments of straight lines between the data, and for computing 
systems without compilers. Other curves will require other 
20 subroutines, but the functions of the subroutines remain the 
same regardless of the curve. 

3. Organization of the Data: Because giving a 
different name to each data point may prove impractical , a 

25 two-dimensional array X may prove effective to store and 

manipulate the data, for instance, with the convention that: 

x i,j :x> <*j>i 

so that the j-th column of the array X contains the point 
30 Xj - (X 1#j ,X 2# j,X 3# j) . 

4. Subroutines for Arithmetic with Vectors: The 
following subroutines may prove efficient for an algorithm 
selected to fit a helix to points in space. 

35 



4.1. SAXPY(a,p,q) . 
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The acronym "SAXPY" stands for "Scalar •A* times X 
Plus Y, H which constitutes a mnemonic device to remember that 
the subroutine SAXPY takes one number a and two vectors x and 
y , and computes (a*x) - y. The vectors x and y may also be 
5 columns of X. In the particular case a = ±1, the subroutine 
gives x±y. Specifically, SAXPY(a,p,q) is again a vector with 
three components, 

SAXPY (a, p,^:- a • p x - q x , 
io SAXPY(a,p7) 2 := a * P 2 - 

SAXPY(a,p,q) 3 : = a • p 3 - q 3 . 

4.2. DOT(p,q). 

The subroutine DOT calculates the dot product of 
15 two vectors: 

DOT(p,q) t-Pxqx + p 2 <fc + P3<*3- 

4.3. SIZE(X) . 

The subroutine "SIZE" computes the Euclidean norm 
20 |xfl 2 , also called length, of a vector x=(x 1 ,x 2 ,x 3 ) in space. 
Specifically, SIZE(x) is a number. 



SIZE(x) := 



4.4. CROSS. 

25 From two vectors, the subroutine CROSS calculates 

their cross product, p q, which has components 
CROSS (p,q)x:=P 2 q 3 - P 3 q 2 
CROSS(p,q) 2 :«p 3 q 1 - p x q 3 
CROSS (P,q) 3 5 == Pi<l2 " P2<Jf 

30 

5. Subroutine DELTA for the Distance Between 
Consecutive Points: The subroutine DELTA (X) calculates the 
Euclidean distance d-j = ||? j+x - Pj| 2 between two points Pj and 
Pj +1 , and the normed difference (Pj +1 - Pj)/dj, which has 
35 length 1, and st res all the r suits into two arrays, wh s 
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nam s may vary, denoted here by D with dimension n and S with 
dimension 3 x n, so that: 

Sj is the j-th column of S, defined by 
Sj:=SAXPY(l f P j+1 , ?j) 

Dj is the length of Sj, 
Dj:=SIZE(Sj) . 

Finally, divide Sj by its length to get a vector of 
10 unit length: 

Sj:«Sj/Dj 

6. Subroutine TNGTS to Compute Tangents: The 
subroutine TNGTS calculates the unit tangent vector at each 
is data point and stores them into two arrays, D of dimension n 
and T of dimension 3 x n: 

T,D:-DELTA(X) • 

Thus, Dj-|3 j+1 -Xj| 2 and Tj«(Xj +1 -j)/Dj for every j € 
20 {l,...,n-l>. 

Other approximations of the tangent vectors with 
finite differences exist. 

25 7. Subroutine CUKVAT for the Normal Vectors and 

Curvatures: The subroutine CUKVAT (T) computes the curvature 
and the normal vector at each data point, and stores the 
results into two arrays, K of dimension n and N of dimension 3 
x n. 

30 N,K:*DELTA(T) 

Thus, Kj-flT j+1 -?j|| 2 and Nj-(T j+1 -?j)/Kj for every j € 
{!,•.. ,n-2>. 

35 8. Subroutine BINORM for the Binormal Vectors: 

The subroutine BINORM calculates the bin rmal vectors and th 
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torsions, and stores the results into two arrays, B of 
dimension 3 x n, and Q of dim nsion n* 

Bj : -GROSS ($^,5.)) , 

Q^DOTCNj.B-,). 

5 

Thus, Bj=TjXNj and Qj - N j • Bj for every j € 
{l,...,n-2). 

9. Subroutine PRMTRS for the Parameters of the 
10 Helix: The subroutine PRMTRS estimates parameters for the 

helix: its curvature KH and torsion QH, and the radius r and 
step s in the parameterization of the helix. 

KH:=(l/[n-2]) (X^-rt^) f 

QH:-(l/[n-2]) (Qi+^Q^) , 
15 c:=KH 2 +QH 2 , 

r:=KH/c, 

s:«QH/c. 

10. Subroutine PLANE for the Plane Perpendicular 
20 to the Axis: The subroutine PLANE fits a plane to the normal 

vectors, for instance, by means of the Total Least Squares 
algorithm described in van Huff el and Vandervalle, supra. 

An alternative consists of using the average 
binomial vector as the direction of the axis, perpendicular to 
25 the plane, stored into an array W with three components: 

W:«(l/[n-2])(Bl + - + B n . 2 ) , 

W:=W/SIZE(W) . 

However, this alternative may prove too inaccurate 
if the data cover only a small portion of one period. Yet 
30 another possibility consists of approximating the normal to 
the plane by the vectors Wj:=Nj x N j+1 for j € {l,...,n-3}, 
and then averaging such vectors: 

W^U/Cn-^lMN! x N 2 + - + N n _ 2 x N n _ 3 ), 

W:=W/SIZE(W) . 

35 Again, the preferred method uses Total Least 

Squares, if it is available. 
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11. Subroutin BASIS f r New Coordinate Axes: The 
subroutine BASIS produces an orthonormal basis (U,V,W) , which 
consists of three mutually perpendicular vectors U,V, and W, 
each with length 1, and such that W coincides with the vector 

5 W from PLANE. 

To this end, consider the unit vector W from PLANE 
as a matrix W with 3 rows and 1 column, and apply to W the QR 
factorization from numerical analysis. The last two columns 
of the matrix Q are the vectors U and V. 

10 

12. Subroutine PRJCTN for the Projections of the 
Data: The subroutine PRJCTN computes the projections 
Pif-rPh of the data onto the plane normal to the axis W, and 
expresses these projections in the system of coordinates with 

15 the basis (U,V,W) . In the plane P spanned by the vectors U 
and V, the projections have coordinates 
Uj-Xj • U = D0T(Xj,U), 

20 which the routine PRJCTN may store in one two-dimensional 
array P, so that P 1#j - u., and P 2#j - Vj. 

13. Subroutine CIRCLE for the Axis of the Helix: 
The subroutine CIRCLE fits a circle to the projections 

25 pl,...,P n in the array P, for instance, with the algorithm of 

Moura & Kitney, supra. 

The circle then determines the axis of the helix, 

which passes through the center (x 0 ,y 0 ,z 0 ) of the fitted 

circle along the direction of W. 
30 The radius of the helix may be either the radius r 0 

of the fitted circle, or the radius r computed above by the 

routine PRMTRS. (For exact fits, the two radii coincide with 

each other.) 

35 14. subroutine ROTATE to Position the Helix: The 

subroutine ROTATE starts a piec f helix at x x with the 
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parameters just described, and focuses on the discrepancy 
between x n and the endpoint of the piece of helix, defined as 
the point where the helix crosses the plane with equation w = 
w . Then the subroutine rotates the helix by an angle <p to 
5 even out the discrepancies at x x and x^ 

15. Summary of the Crude Algorithm: Read the data 
into the array X. 

Call TNGTS. 
10 Call CURVAT. 

Call BINORM. 

Call PRMTRS. 

Call PLANE. 

Call BASIS. 
15 Call PRJCTN. 

Call CIRCLE. 

Call ROTATE. 

Results: An orthonormal basis (U,V,W), with W along 
the axis of the helix, a point <u 0 ,v 0 ,w 0 ) , center of the 
20 fitted circle, where the axis passes, the radius r of the 
helix, the step s of the helix, the angle of rotation v to 
rotate the helix to fit the data better, the constant 
a = l/^r 2 +s a for computational convenience, the curvature KH, 
and the torsion QH, of the helix, which now has the formula: 
25 u(t)»u 0 +r*cos(a-t) , 

v(t)=v 0 +r*sin(a-t) , 

w(t)=sa-t, 

for w^/ts-a) st*w n /(s*a) . 

30 EXAMPLES 
Example 1 

This example demonstrates an embodiment of the 
present invention which is further developed and refined in 
subsequent examples. Stereotactic neurosurgical implantation 
35 was performed with polysulfone (PS) and polyacrylonitrile 

polyvinylchl ride (PAN-PVC) copolymer tubing through most f 
the length of the n ©striatum by the methods of the present 



WO 96/03087 



PCT/US95/09391' 



99 

inv nti n. The planned trajectory was confirm d by x-ray. A 
small amount of dy was allow d to diffuse into the brain all 
along the course of the tubing, confirming the circumfusion of 
the dye and the positioning of the tubing by examination of 
histological sections. 

Post-mortem cranial and brain specimens from two 
adult female pigtailed macaques (Macaca newegtrina) were 
acquired from the Tissue Distribution Program of the 
University of Washington Regional Primate Research Center. 
Euthanasia was performed by Nembutal overdose in the manner 
approved by the university. 

A spreadsheet was devised to provide specifications 
for passing the needle through a given brain structure. 
Inputs to the spreadsheet included (1) the stereotactic 
coordinates (AP, ML, and DV) of three points to be included in 
the trajectory, (2) the guide-tube length and inner diameter, 
and (3) the outer diameter of the curved needle. The 
spreadsheet also provided a choice of species (longtailed 
macaque, pigtailed macaque, or baboon) by incorporating the 
locations of the anterior commissure in each species as a 
point of reference, as pictured in standard brain atlases. 
Formulations for adjusting the stereotactic position on the 
basis of ventriculography were also incorporated in the 
spreadsheet, but were not used for the post-mortem test, for 
which ventriculography would be impractical. . 

Based upon these inputs, the spreadsheet provided 
the diameter of curvature, the angle and position at which the 
guide-tube was to be inserted, the points at which the tubing 
would intersect coronal planes, and the length of curved 
tubing that would reach from the top of the guide-tube to the 
point at which it re-emerged from the brain. The spreadsheet 
was written in Excel for Windows in the Windows 3.0 
environment (Microsoft Corp., Redmond WA) , implemented on a 
Dell System 310 IBM-compatible 80386 microcomputer (Dell 
Computer Corp, Austin TX) • 

Th spreadsheet was cued by entry of th length (35 
mm) and int rnal diam ter (0.965 mm) of the guid -tube 
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(18 -gauge thinvall) , the outer diameter (0.330 mm) of the 
curv d needle (29-gauge) , th choice of species and reference 
atlas, and the stereotactic coordinates of three non-collinear 
points on the planned trajectory. From these coordinates, a 
5 series of automatic matrix operations in the spreadsheet 

derived the orientation of the plane (plane "PP") containing 
the three points to the standard stereotactic coordinates, and 
the coordinates of the center of the circle in PP specified by 
those three points. The center was determined by taking the 

10 three original points as the vertices of a triangle, for which 
the perpendicular bisector planes of each side must intersect 
in a line perpendicular to the plane PP and equidistant from 
the three points. The intersection of the plane PP with this 
line (i.e., with any two of these bisector planes) was the 

15 center of the circle in PP determined by the points. 

The two intersection points of the circle and any 
given coronal plane were determined from the equations for PP, 
the circle, and the coronal plane. (For the planned 
trajectory, with the straight guide-tube lying in the most 

20 posterior coronal plane, the upper point was in the brain only 
for the most anterior coronal planes • ) The direction cosines 
of the tangent (in PP) to the circle, at each of these 
intersection points, were determined from the intersection of 
PP with the plane perpendicular to PP and tangent to the 

25 circle at the given point. For a guide-tube of the given 
length placed at the angle specified by these direction 
cosines, the coordinates of the top and the tip were derived 
from the direction cosines, taking into account the fact that 
the curvature of the needle would begin inside the guide-tube 

30 at a calculable distance above its tip. These calculations 

were carried out for all coronal planes represented in the atlas. 

The most posterior point on the circle was 
calculated numerically by substituting values to determine the 
minimal AP level at which there were real ML and DV values. 

35 At this point, the angle of the tangent-plane with the AP axis 
in th sagittal plane must b 90°. Th spreadsheet provided 
the coordinates of the top of the guide-tube when th tip was 
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inserted to this point, allowing confirmation that the guide- 
tube length was suff ici nt to clear the top of the skull and 
allow for insertion of the circular needle. 

The spreadsheet also provided two other coronal 
levels, at which the trajectory was expected to lie (1) 
outside the brain for the first tine and (2) furthest outside 
the brain, when the needle was maximally inserted in the 
guide-tube. The spreadsheet provided the length of tubing 
required for the trajectory, including the portion above and 
within the guide-tube, as well as the portion circling through 
the brain and extending outside the brain at maximal 
insertion. 

The putamen was targeted for implantation by the 
preparations. The head was positioned in the stereotactic 
instrument and scalp and bone tissue were removed from the 
working area. The carrier was moved to the anteroposterior 
(AP) position over the post-central gyrus specified by the 
spreadsheet, the mediolateral (ML) scale adjustment was made 
as specified, and the guide-tube, with a stylet inserted, was 
lowered by the stereotactic dorsoventral (DV) adjustment to 
the specified DV scale position (final location of tip: P2, 
Rll.l, D29.7, relative to earbar zero). The stylet was 
removed and the curved needle inserted. A 2 9 -gauge needle, 
curved to the specified radius, was lowered through the guide- 
tube, and as it passed out the tip of the guide-tube and 
through the brain, the plane of its emerging arc was kept 
parallel to the side of the rectangular needle-holder. When 
its tip emerged through the surface of the prefrontal cortex 
(about A25) , it was inserted further into the guide-tube to 
bring the tip about 3 mm out of the brain. PS tubing (Amicon, 
Danvers MA; MW cutoff 100,000 Da, OD 0.5 mm, from H1P100-20 
Diaf lo Hollow Fiber Cartridge) , or PAH-PVC tubing (MW cutoff 
50,000 Da, OD 1.0 mm), previously cleared of glycerine by 
ultrafiltration with physiological saline and kept moist and 
filled in a saline bath, was fastened to the tip of the needle 
with Vetbond (Tissu Adhesive, N. 1469; 3M Animal Care 
Products, St. Paul MN) . The needl was slowly withdrawn, 
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r tracing its course, pulling the polymer tubing into 
prefrontal cortex, down through the length of the neostriatum, 
and up to the lower tip of the guide-tube. The guide- tube, 
containing the curved needle and pulling the polymer tubing at 
its tip, was then raised by the stereotactic DV adjustment. 

The radius of curvature of the needle, the angle at 
which the guide-tube was inserted, the points at which the 
tubing intersected coronal planes, and other useful parameters 
were provided by customized spreadsheet software as described 
above. 

The 29-gauge insertion needle was curved to the 
diameter specified by the spreadsheet by (1) wrapping needle 
stock several times in a tight helix around a mandrill rod of 
about 0.25 the desired diameter, (2) passing the resultant 
coil through the guide-tube, and (3) manually adjusting the 
helical displacement to approach a circular contour. 

A Kopf stereotactic instrument (#1404, David Kopf 
Instruments, Tujunga CA) was set with the rotary scale on the 
carrier facing squarely anterior, so that at a given AP level, 
the stereotactic shaft (mounted by the carrier to the 
stereotactic arm), along which DV adjustment is normally made, 
could move only in a plane vertical to the stereotactic 
horizontal. Because a guide-tube oriented tangential to the 
most posterior point of the circle must have a vertical trace, 
as viewed from the side, this placement pre-determined the 
first of three components of the proper setting for the plane 
PP. 

A custom-designed needle-holder, in the form of a 
rectangular stainless-steel block, was mounted to the lower 
end of the stereotactic shaft by a universal ball- joint. The 
ball-joint was first adjusted to bring the needle parallel to 
the stereotactic shaft along which DV adjustment is normally 
made; this was accomplished by setting this shaft to vertical, 
loosening the ball joint, placing a rectangular, precision- 
machined, stainless-steel calibration bar squarely across the 
arms of the stereotactic instrument, adjusting the DV s tting 
until th needl -h Ider sat squarely on the bar, then 
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tightening the ball joint. The ster otactic shaft was then 
rotated by an angle prescribed by the spreadsheet and drawn 
with a protractor on the calibration bar; the plane PP now 
stood orthogonal to the stereotactic horizontal but at an 
5 acute angle to the stereotactic sagittal plane. With the 

shaft tightened in place, the third component was provided by 
adjusting the rotary scale on the carrier to an angle 
prescribed by the spreadsheet. With the angles fully 
adjusted, the tip of the guide-tube in the needle holder was 

10 brought to earbar zero, and the AP, ML, DV settings were 
recorded in the spreadsheet. 

A ventriculography could have been performed in a 
living monkey, to bring the intercommissural line parallel to 
stereotactic horizontal and to locate the anterior commissure 

15 within the adjusted coordinate system with greater accuracy 

than available from the atlas alone. The spreadsheet allowed 
the necessary inputs to be made, but for the post-mortem test 
material, this step was omitted and the atlas values were 
considered the best estimates available. 

20 The spreadsheet combined the AP, ML, and DV earbar 

zero scale readings, the estimate of the position of the 
anterior commissure, and the location of the point most 
posterior on the circle, to provide the AP, ML, and DV 
settings for lowering the guide-tube to the proper position. 

25 With both ends protruding from the brain, the 

tubing was filled with radio-opaque iohexol (Sterling- 
Winthrop) and a lateral x-ray was taken to confirm the 
position and contour of the implant. The iohexol was then 
replaced with a 0.75% solution of malachite green dye. The 

30 dye was replaced with fresh dye every 30 min for 1.5 h (PAN- 

PVC tubing) or 3 h (PS tubing), to test for diffusion into the 
brain parenchyma surrounding the tubing. Finally, lot 
formalin was introduced into the tubing, the brain was removed 
and blocked, and the portion containing the tubing was fixed 

35 by immersion in formalin for 24 h, then cryoprotected by 

imm rsion in 30% sucrose for 3 days. Coronal froz n sections 
wer cut at 1-mm intervals through ut th course of the 
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implant and mounted on glass slides. One series was compared 
with th 1-mm intervals in the brain atlas and the best match 
for each section was selected. The position of the implant in 
the section was transferred by reference to the comparable 
5 neighboring structures in the atlas, and the ML and DV levels 
of that position were then read from the atlas. The measured 
readings for each section were compared with the planned 
position for that section provided by the spreadsheet. 

The implanted tubing was located in the dorsal half 

10 of the neostriatum for approximately 16 mm, reaching 80% of 
the AP extent from the anterior tip of the caudate to the 
posterior tip of the putamen. A halo of green dye surrounded 
the tubing in each section, fading at the edges, with a 
diameter of about 2 mm. 

15 X-rays displayed a smooth course with a shape very 

similar to that predicted for a lateral view of the planned 
trajectory. The sequence of 1-mm coronal sections matched 
well with sections in the atlas. Within the measured series 
both the posterior and anterior limits of the trajectory were 

20 located about 2 mm posterior to the expected positions 
relative to the anterior commissure. For the other two 
dimensions, measurements from histological slides (PS tubing) 
indicated a trajectory generally more dorsal by about 2 mm and 
slightly more medial than planned at many points. 

25 This study has demonstrated that permselective 

polymeric tubing can be accurately implanted over a 
substantial distance in a primate brain, leaving both ends 
available for the introduction .of biochemicals, which are then 
free to diffuse out of the tubing and into the brain 

30 parenchyma. In tissue from the test subjects, the implant 
passed through most of the length of the neostriatum. The 
x-ray indicated a smooth trajectory for the polymer tubing, as 
planned, suggesting that irregularities in the rough numerical 
ML and DV estimates, made from the mounted sections, were due 

35 largely to measurement error associated with routine 

histological procedures. The measurements from histological 
slides indicated non-random systematic err rs in the 1 cation 
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of the trajectory, attributable to the variation of the 
individual animal from the estimate provided by th atlas. 
Placement was reasonably accurate, but the relative positions 
of bony and brain landmarks can vary substantially, making it 
necessary in nonhuman primates to improve stereotactic 
accuracy by performing ventriculography or magnetic resonance 
imaging in the living animal. Such a procedure was not 
possible for this post-mortem tissue, so that small errors in 
the placement of the circular trajectory were to be expected. 
The apparent deviations suggest that in this tissue, as 
positioned in the stereotactic instrument, brain landmarks 
were located 2 mm anterior, 2 mm ventral, and 0.8 mm to the 
right of their atlas locations. Compensation for these errors 
due to the deviation of this subject from the atlas series 
would have further improved the accuracy of placement. 

The green haloes around the tubing in tissue 
sections were of dimensions comparable to those seen after 
point injections of malachite green in rodents (Myers, 
physiol. Behav. . 1, p. 171, 1966). This suggests that passive 
diffusion is ample for passage of a biochemical, highly 
concentrated in the tubing, into surrounding brain tissue, 
even post-mortem. Beyond this observation, the distance 
covered by this diffusion over 1.5-3 h, in post-mortem tissue, 
with occasional replacement of the fluid in the tubing from a 
reservoir, cannot be assumed to represent the distance that 
would be covered by an endogenous biochemical such as dopamine 
applied chronically in a similar manner. This distance should 
be expected to depend not only on the concentration, the 
number of hours or days of treatment, and the rate of 
replacement, but also on the presence of co- infused metabolic 
inhibitors, the geometry and fluid dynamics of the 
extracellular space in gray and white matter, and long-term 
brain biochemical adaptations that may occur in response to 
the chronic presence of a source of biochemical. 

PAN-PVC has been shown to be highly biocompatible 
with brain tissu , with inf lammati n subsiding aft r a few 
weeks and th astrocytic zon around the implant d cr asing 
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from over 400 fim to under 40 /an over 12 we ks (Winn et al., 
Exoer. Neurol. . 105, pp. 244-250, 1989). The chronic presence 
of a stable high concentration of a biochemical in a 
renewable, biocompatible source that is well -integrated into 
5 the brain should lead to a steady state in which much of the 
neostriatum can be profoundly influenced from a single 
implant. In rats with unilateral substantia nigra lesions, 
infusion of dopamine by implanted cannula over several days 
into brain tissue has been shown to have substantial effects 

10 on apomorphine- induced rotational behavior and on striatal 
levels of dopamine, while intraventricular infusion at the 
same rate had no behavioral effect and raised dopamine levels 
only within 1 mm of the ventricle (Kroin et al., J. 
Neurosurcr. . 74, pp. 105-111, 1991). Chronic infusion of 

15 labeled dopamine by a standard cannula over several days in 
rabbits has led to measurable label as far as 10 mm from the 
cannula tip (Sendelbeck and Urguhart, Brain Res. . 328, pp. 
251-258, 1985). In an earlier study (Dubach, Neuroscience . 
45, pp. 103-115, 1991) it was shown that acute point 

20 injections of highly concentrated dopamine, along with an 

inhibitor of monoamine oxidase, can leave high tissue levels 
elevated for as long as two hours after injection, over a 
region that extends more than 3 mm from the center of the 
site. 

25 

Example 2 

Striatal implantations and subsequent rotational 
observations were made in two baboons ( Papio cvnocephalus ) , 
and the parahippocampal region has also been implanted in a 

30 long-tailed macaque ( Macaca fascicular is ) • 

The baboons were treated unilaterally with the 
dopamine neurotoxin MFTP via the carotid artery several months 
before implantation using the technique originated by 
Dankiewicz et al. ( Life Sci. . 39, pp. 7-16, 1986). Flow was 

35 directed from the common carotid up into the internal carotid 
during th infusion by occluding the external carotid. The 
internal carotid was easily distinguishable by its deep 
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passage into the neck in the bab on; the two more superficial 
branches of the external carotid artery typically bifurcated 
immediately upon leaving the common carotid. The infusion 
rate (60 cc over 12 minutes) was slightly faster than that of 
Bankiewicz (60 cc over 15 minutes) (Bankievicz et al., supra). 
Isoflurane anesthesia (Kor dower et al., J- Neurosura. . 73, pp. 
418-428, 1990) was used for the lesion surgery. A high 
concentration of MPTP (1 mg/kg, 2.5 times the concentration 
used by Bankiewicz) was required for consistent results in 
baboons; possibly this is because the brain/body ratio is 
larger, so that to get the same amount to the brain requires a 
larger "per kilogram body weight" amount. 

Calculations for planning trajectories were 
organized in the form of a spreadsheet implemented under 
Microsoft Excel (Microsoft Corp* , Redmond WA) . The 
spreadsheet was augmented by access to magnetic resonance 
image (MRI: General Electric 1.5 Tesla) files manipulated by 
Microsoft Visual Basic and displayed by Adobe Photoshop in a 
Microsoft Windows environment. Inputs included the dimensions 
of the tubing and wire, landmarks that established the 
orientation of the brain in a coordinate system shared by the 
MRI and the stereotactic instrument, and three points through 
which the implant was intended to pass. The spreadsheet 
derived the orientation of the plane and the center and radius 
of the circle in that plane that were specified by the three 
"target" points. From these results it determined the 
intersection points of the circle with each coronal plane, if 
any, and the direction cosines of the tangent to the circle at 
each intersection point. Two tangent locations on the circle 
(for the guide-tube cannula and for the emergent straight 
wire) were selected by the user by choosing two of these 
coronal planes* Given the tangent selected as the course of 
the guide-tube cannula, the spreadsheet determined the 
coordinates of the tip of the cannula at full insertion 
(taking into account the curvature of the tubing that would 
remain ins id the tip of th cannula) s that the us r could 
zer the st reotactic equipment appr priat iy. Th user 
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entered the zero s ttings in the spreadsheet. The user also 
specified th coronal plan s at which the trajectory would be 
expected to enter and exit the target structure, and the 
coronal planes at which the tangents (cannula and wire) would 
5 be expected to intersect the cranium. Output cells in the 

spreadsheet then reflected the stereotactic angle settings and 
Vernier-scale settings for the guide-tube cannula, and the 
required lengths of nitinol tubing, nitinol wire, and 
permselective implant tubing. 

io MR images were made for each monkey as an aid to 

stereotactic placement. A non-ferrous stereotactic instrument 
was designed and built for this and other projects, in the 
Regional Primate Research Center Machine Shop. The "spoiled 
GRASS" format on the GE 1.5-tesla magnet was used for making 

15 these images. This format provides a "voxel " (3-D pixel) that 
is 0.8 mm on a side, resulting in an accurate, continuous, 
customized "atlas" of each monkey's brain. 

The spreadsheet provided several options for 
planning various types of trajectories, including a Graphical 

20 User Interface option, by which as many as 30 coronal MR 

images were displayed simultaneously on the monitor screen. 
Individual images could be enlarged and the location of 
landmark and target points entered graphically by a mouse. 
The system transferred these coordinates to the spreadsheet, 

25 which calculated the trajectory passing through the target 
points, and projected onto each of the 30 coronal views a 
marker for the point at which the tubing would pass through 
that section. This set of views greatly simplified evaluation 
and adjustment of a putative trajectory and determination of 

30 the planes near which that trajectory entered and exited the 
cranium and the target structure. 

The putamen (or the parahippocampal cortex in the 
longtailed macaque) was targeted for implantation by computer 
methods outlined above. The implantation consisted of 

35 establishing a curvilinear course down into and back up out of 
th brain with a shaped m tal 1 ad r, and th n drawing a piece 
of flexible tubing in th revers dir ction back thr ugh the 
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brain along that course. A t 1 scoping set of concentric 
lengths of metal tubing and wir traveled al ng the planned 
course down into the brain. A stainless-steel guide-tube 
cannula (21-gauge encased in 18 -gauge) passed through a 
5 burr-hole in the cranium and entered the brain at a prescribed 
point on the cortical surface and followed a straight line 
down to a "tangent point. H At that point a radially curved 
span of nitinol tubing (approximately 2 6 -gauge) emerged to 
trace an arc in the depths of the brain around to another 

10 prescribed point. At that point a length of fine nitinol wire 
(10-mil) emerged tangent i ally from the curved tubing and 
traveled up out of the brain along a straight line, to pierce 
through the cortical surface from underneath at a point not 
far from the original entry point. Nitinol, a "superelastic" 

15 nickel -titanium alloy, permits a much tighter curvature than 
stainless steel, as well as enabling the use of the fine wire 
to trace the second straight leg of a "hairpin" course through 
the brain. The cannula, curved nitinol tubing, and straight 
nitinol wire were mounted by set screws respectively into 

20 lower, middle, and upper stainless-steel collar mounts 

(disks) , which were kept in register by similarly mounted 
guide rails. The guide-tube cannula extending from the lowest 
disk was mounted by set screw into a stainless-steel block 
fastened to a needle carrier on a conventional Kopf 

25 stereotactic instrument (#1404, David Kopf Instruments, 

Tujunga CA) . During insertion, the final passage of the wire 
was stopped short of the underside of the cranium, x-rays were 
taken to confirm the point at which it would emerge, and 
another burr-hole was made in the cranium. 

30 When the metallic elements were fully inserted, 

with the nitinol wire emerging from the brain, a length of 
polyethylene tubing (PE-10) was fastened to the wire with 
cyanoacrylic glue (Vetbond: 'Tissue Adhesive, N. 1469; 3M 
Animal Care Products, St. Paul MN) and the disks anchoring the 

35 wire and tubes were raised in series, pulling the PE-10 down 
into and back up out of the brain, f 11 owing the same 
traject ry in reverse rd r. The PE-10 tubing had previously 
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been glu d into a length of permselective polyacrylonitrile 
polyvinylchloride (PAN-PVC) copolymer (MW cutoff 50,000 Da, OD 
1.0 nun) tubing, in such a way as to provide a patent span of 
PAN-PVC in the targeted region of the brain and non-patent 
segments lined with PE-10 extending from this span up out of 
the brain; the span lengths were all read from the 
spreadsheet. The construction of the compound PE/PAN-PVC 
tubing was accomplished by (1) first slipping a length of PAN- 
PVC (one of the occluded segments) all the way onto the PE-10, 
(2) then inserting the tip of the PE-10 just inside another 
length of PAN-PVC (the patent span), (3) placing a small 
amount of Vetbond on a mark 5 mm up the PE-10, and (4) pushing 
it into the PAN-PVC tubing up to the mark, (5) then sliding 
the previously mounted PAN-PVC segment down to the glue joint. 
Another length of PE-10 was connected with the patent span of 
PAN-PVC at the other and in the same manner. This compound 
flexible tubing was then pulled far enough into the brain that 
the patent semi -permeable portion was exposed to the target 
area and the internally reinforced (PE-10) and externally 
biocompatible (PAN-PVC) tubing passed through remaining brain 
areas and the cranium. Tubing was bathed in saline as it 
passed into the brain. The friction at the brain/ tubing 
interface was small enough, and the support provided by the 
brain as the tubing passed around the radial curve was large 
enough, that the PAN-PVC/PE-10 tubing traveled along the 
course established by the nitinol and stainless steel, which 
pierced but not did not tear the brain tissue. 

The burr-holes in the cranium, which provided 
access to the brain for the tubing, were repaired by placing a 
thin layer of gel foam at the depth, applying Vetbond to the 
sides of the hole and to 1-2 mm of the cranial surface 
surrounding the hole, and then introducing a few drops of 
Kooliner acrylic. The PAN-PVC tubing (reinforced with PE-10) 
that emerged from the brain was thus joined firmly to the 
cranium. A length of PE-50 tubing, previously sleeved in 
silic n tubing for tissue c mpatibility and passed 
subcutaneously from a point f emerg nee on the back up t 
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cranial incisions, was joined with Vetbond and silk sutur to 
the PE-10. This continuous length of tubing was then sutured 
to the scalp during closure, so that the entire 
implant — acrylic and compound PAN-FVC/PE tubing — was 
subcutaneous down to the point of emergence on the back. 

The monkeys were equipped with nylon jackets and 
flexible-cable tethers for drug administration. The two 
lengths of PE-50 encased in silicone tubing, from the two ends 
of the implant, provided for safe subcutaneous passage and 
also continued up the flexible tubing fastened to the monkey's 
nylon jacket and then into the rotating shaft of the rotometer 
above the cage, where they were led out through a hole in the 
shaft and connected via three-way stopcocks to a pair of 1-cc 
syringes fastened to the shaft to rotate along with it and 
along with the monkey below. The stopcocks permitted repeated 
emptying and filling of the syringes without disturbing the 
continuity of fluid in the access tubing that passed through 
the brain. The "outflow" syringe, connected to the 
posterior-cranial PE-50 tubing, was mounted upside down; 
during replacement of fluid in the access tubing, its plunger 
was interlocked with that of the "inflow" syringe just below 
on the shaft, to provide a balanced push-pull action and to 
avoid any unnecessary pressurizing or decompression of the 
fluid in the access tubing in the brain. 

Adaptations in the tether system enabled it to 
serve as the drive-shaft for baboon rotometry as well as the 
conduit for drug administration. The cable was strengthened 
and stiffened by a wire-braid wrapping. The connection to the 
jacket for Monkey 2 was also set at a right angle and 
reinforced with a sturdy stainless-steel housing that 
prevented its destruction, after the loss of several earlier 
implementations due to the challenging motor behavior of this 
monkey. A channel in the wire mesh of the top of the cage, 
from the middle to the front, allowed for transferring the 
monkey, tether cable, and syringe/rotometer shaft as a single 
unit from one cage t anoth r, t permit removing and cleaning 
the cages. 
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Rotational behavior was monitored continuously by 
rotometer methods, and the effects of lesions and implants 
were studied without the peripheral administration of CNS 
stimulant drugs. The behavior of lesioned animals normally 
5 involves some contralateral turning as well as ipsilateral, 
but the distribution of this "paradoxical" turning was not 
evenly spaced through the day. In both baboons, the greatest 
amount of paradoxical turning occurred during an active period 
in the middle of the day, for the 2-3 hours preceding feeding. 

10 Drug applications via the access tubing were therefore made in 
the morning, before this period. 

A concentrated solution of apomorphine (0.0045 M) 
was placed in the access tubing of Monkey 1 on two occasions, 
separated by 11 days. The first application consisted of 

15 filling the entire length of tubing (approximately 7 feet) 
with a concentrated apomorphine solution by the push-pull 
action of the syringes described above, then obstructing the 
outflow line and slowly pushing in 0.25 cc under pressure. 
The only patency in the tubing was the 19-mm permselective 

20 portion in the putamen, the brain structure most affected by 
the MPTP lesion and by Parkinson's disease in humans. The 
second application was more gentle, by diffusion alone, using 
the same high concentration of apomorphine, but without any 
pressurized injection. A break in the line in the flexible 

25 cable a few days after the second treatment prevented further 
testing of Monkey 1. Treatments in Monkey 2 began with a low 
dose of apomorphine 1/100 of the original dose. Fluid at this 
concentration was placed in the implant on the second day and 
again on the third day after the surgery. A higher dose (1/50 

30 original) was placed on each of the next three days, then 

1/20, 1/10, and 1/3 for two days each, and finally the full 
concentration for the final day, after which saline was 
returned to the line. 

Two days later, on the premise that oxidized 

35 apomorphine might be occluding the molecular pores in the 

tubing, the pr ssure resisting infusion of fluid was tested by 
releasing the syringe n edle n th upper lead from its 
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syringe and stopcock, allowing fluid to infuse under the 
pr ssure of gravity, and by measuring th rate at which the 
fluid level dropped in the tubing, over a distance of 3 mm. 
By comparing this to the rate of fall under similar conditions 
in vitro with the patent span of tubing resting in saline at 
atmospheric pressure, an estimate of "circumtubular* pressure 
was calculated. Subsequently, an attempt was made to clear 
potential occlusion of the pores with a solubilizing agent, 
dimethyl sulfoxide (DHSO) . DMSO (10% in normal saline) was 
introduced and left in the tubing for 30 minutes, then 
replaced with saline. The drop-rate was then re-measured, 
then DMSO was introduced again for another 60 minutes, and the 
drop-rate was measured a third time. 

Nine days after the last apomorphine exposure, 
dopamine and tranylcypromine were presented together in 
solution and left in place without replacement for the next 
four days. At this point a break in one line outside the 
brain required that behavioral experiments be terminated. 

The tubing implanted in Monkey 2 was visualized by 
MR imaging after behavioral experiments were terminated. 
Series of images were made just before and just after the 
application of gadolinium (3 parts in 50 parts saline) , as 
well as 25 min and 50 min after gadolinium. The tubing 
implanted in Monkey 1 was also visualized, but without 
gadolinium. 

The portion of the brain containing the tubing was 
fixed by immersion in formalin for 24 h, then cryoprotected by 
immersion in 30% sucrose for 3 days. Coronal frozen sections 
were cut throughout the course of the implant , and Nissl 
sections and immunoh 1st ochemica 1 sections (tyrosine 
hydroxylase and glial fibrillary acidic protein) were prepared 
at approximately 1-mm intervals by methods previously 
described. 

The pressurized infusion of 250 fil of concentrated 
apomorphine had both acute and lasting effects in Monkey 1. 
Th m nk y almost imm diat ly (2 min) began turning in th 
direct! n contralat ral to the side f the 1 si n and the 
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implant. This turning was unlik the typical contralat ral 
turning for this monkey. Typical contralateral turning was 
performed bipedal ly or on all fours with the face near the 
walls of the cage, directed straight ahead or toward the 
ipsilateral side. The "topography" of the apomorphine-induced 
turning behavior, on the other hand, was comparable with that 
of turning induced in lesioned rats by peripherally 
administered apomorphine. The face was usually directed 
toward the contralateral shoulder. Turns were fast and 
"tight" and were made with minimal use of the hind limbs, 
especially the contralateral hindlimb, which was used as a 
pivot at an early point, within 3 min of drug infusion. Later 
on, the hindlimbs more often were flexed or collapsed. Some 
of the behavior could be described as "ischial spinning", 
mixed with bursts of bipedal spins or jumps (within 5 min) . 
Turning was less continual by this time, and most turns were 
single full rotations, explosive movements at intervals of 
5-40 seconds, with a prominent weaving and crouching of the 
upper body between turns. These effects were so intense that 
the apomorphine was replaced with saline (17 minutes after 
treatment) . During this procedure bouts of multiple ischial 
spins occurred (21 minutes after treatment) , and such bouts 
and isolated rapid turns continued to occur for about an hour, 
gradually becoming slower and more often aborted to a 
head-and-shoulder gesture in that period of time. Fig. 11A 
(apomorphine under pressure) represents contralateral turns 
per minute for each minute through the day, with the time 
marked on the x-axis in hours. 

in a second trial after 11 days to assure recovery, 
introduction of apomorphine at the same high concentration, 
but without pressurized infusion, also produced ischial 
spinning of a similar topography and latency of onset. The 
rotation was somewhat less "explosive" in character, but more 
continual; it involved less ataxic weaving of the upper body, 
but more scratching, which made use of the affected 
c ntralat ral limb, which was n t normally used in this 
manner. Again, ffects app ar d so rapidly and intens ly that 
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an attest was made to terminate them by r placing apomorphine 
in the tubing with saline (18 min after tr atment) . In 
addition, 28 minutes later, the saline was replaced by the 
dopamine antagonist droperidol (0.0053 M) . immediately after 
this exchange the monkey became very inactive. Although the 
droperidol was left in the tubing for 7 hours, however, it did 
not prevent the midday bout of contralateral turning that was 
typical for this monkey. Pig. 11B (apomorphine by diffusion) 
shows the contralateral turns in response to the treatment. 

Monkey 2 was implanted twice. The stabilization of 
the tubing as it emerged from the cranium after the first 
implantation was inadequate because Vetbond was omitted during 
preparation of the posterior acrylic pad, which quickly broke 
loose, interrupting the line. Examination of the preparation 
and removal of the implant after this occurrence, under 
inevitably non-sterile conditions, may have led to bacterial 
contamination of cranial and brain tissue, which was treated 
prophylactically with antibiotic. Six weeks later, a second 
implant was introduced, with the pads stabilized with 
additional Vetbond between the cranium and the acrylic pad as 
described above. This implant remained serviceable for five 
weeks. 

For a 12-day period beginning two days after the 
second implant surgery, the fluid in the tubing was changed 
every day, gradually increasing the dose of apomorphine from 
1% to 100% of the previous 0.0045 M concentration. Prom a 
starting point of 5-10% contralateral turning per day, the 
monkey reached an average of 21% contralateral turning over 
the last week of this period, under the influence of the 
higher doses, with a peak day at 35%. Replacement of the 
highest concentration with saline alone, after this period, 
resulted in a gradual decline back to baseline over 7-8 days. 

During this period the pressure resisting an 
infusion of fluid was found to be well beyond the normal 
limits of intracranial pressure. Given the possibility that 
the mol cular por s in the tubing had been larg ly occluded by 
chronic exposure t oxidized apomorphine, an attempt was mad 
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to clear any such occlusion. The solubilizing agent DMSO was 
administered, first for 30 minut s, and then for 60 minutes, 
and the drop-rate increased after each exposure. 

After a total of nine days at about 10% 
contralateral turning, including the day of these pressure 
tests, a solution of dopamine itself (0.045 M) was placed in 
the implanted tubing, along with tranylcypromine (0.0188 M) , a 
monoamine oxidase inhibitor to delay the metabolic disposal of 
the dopamine. The day after this change, the monkey showed 
31% contralateral turning, but returned on the following day 
to baseline levels. The fluid was not exchanged and the 
intactness of the line was not examined until four days after 
placement of the dopamine/ tranylcypromine solution, at which 
time the outflow line was found to be interrupted. 

Histological examination of Monkey 1 indicated the 
placement of the tubing in the medial putamen, internal 
capsule, and anterior caudate nucleus throughout most of the 
length of the neostriatum, and indicated that the lesion had 
greatly reduced TH immunoreactivity in the striatal dopamine 
terminal areas. 6FAP immunoreactivity indicated activation of 
glial cell somata in the vicinity of the tubing, where labeled 
cells were approximately three times more common than in areas 
distant from the tubing. The trabecular walls of the tubing 
which surround the inner permselective membrane were marked at 
some points by intrusions of labeled material. Control 
sections reacted with normal rabbit serum as the primary, or 
with vehicle alone without any primary antibody, displayed 
different patterns of peroxidation in the walls of the tubing. 
It was not possible by IM observation alone to determine to 
what extent glial extensions had invested the walls of the 
tubing, or to what extent the "labeled" objects consisted of 
oxidized apomorphine. It is interesting to note, however, 
that every section subjected to free-floating 
immunohistochemical processing consisting of more than 24 
hours of shaking still retained its PAN-PVC tubing 
cross-s cti n. 
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In vivo MRI evaluati n f th tubing in M nk y 2 
indicated that its course through the brain clos ly matched 
the planned trajectory. The tubing when filled with saline 
appeared as a few dark voxels, an absence of signal similar to 
5 the lateral ventricle. The introduction of "caged" gadolinium 
into the tubing resulted in a strong signal within the tubing, 
providing a detailed assessment of its trajectory. The second 
and third MRI views at 25 and 50 minutes after introduction of 
gadolinium indicated a gradual spread of signal from the body 

10 of the tubing into the surrounding tissue. This could be 

observed by close comparison of the original pictures, but was 
more clearly evident from subtractions performed within the 
Adobe Photoshop environment. 

Histological examination of Monkey 2 indicated the 

15 presence of two tubular implants extending through most of the 
length of the put amen. It was possible to identify the later 
implant in the more lateral of these two positions, 
corresponding to the gadolinium-marked tubing of the In vivo 
MRI evaluation, and the earlier implant, which was never used, 

20 in the position 2-3 mm more medial. The later implant, which 
was used chronically for applying dopamine and apomorphine, 
and acutely for applying DMSO and gadolinium, was associated 
with evidence of an intense local inflammatory response, but 
the medial implant did not appear to be directly associated 

25 with this. A reasonable interpretation would be that the 

presence of the tubular implant itself was innocuous, but that 
one or more of the substances applied via the tubing is 
capable of causing inflammation in the concentrations 
employed. 

30 In vivo evaluation of the parahippocampal implant 

in Monkey 3 by MRI indicated that the tubing followed a course 
very similar to that planned. The trajectory followed a 
straight and nearly vertical path of entry posteriorly, as 
planned, down to a point just behind the anterior enlargement 

35 of the inferior horn of the third ventricle, then passed under 
the ventricl and hipp campus, in the parahipp campal area, to 
an anterior point just medial t the ventricle, fr m which it 
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ffflftT T fi CTATMED IS; 
X 1. A devic having a proximal end and a distal 

2 end for percutaneous localization of a foreign object in a 

3 body, comprising: 

4 a plate-like base collar mount having a proximal 

5 surface, a distal surface rigidly affixed to at least one 

6 perpendicular guide rail, and a channel between the upper and 

7 lower surfaces through which channel a guide-tube cannula 

8 extends from the proximal surface of the base collar; 

9 a plate-like first movable collar mount slidably 

10 located on each guide rail rigidly affixed to a curved first 

11 extendable cannula having a proximal end and a distal end, 

12 wherein the distal end of the first extendable cannula is 

13 affixed to the first movable collar in a channel through the 

14 first movable collar and the proximal portion of the first 

15 extendable cannula is in the guide-tube cannula; and 

16 a plate-like second movable collar mount slidably 

17 located on each guide rail rigidly affixed to a second 

18 extendable cannula or wire having a proximal end and a distal 

19 end, wherein the distal end of the second extendable cannula 

20 or wire is affixed to the second movable collar in a channel 

21 through the second movable collar and the proximal portion of 

22 the second extendable cannula or wire is in the first 

23 extendable cannula; 

24 wherein proximal sliding of the first movable 

25 collar will extend the first extendable cannula through the 

26 proximal end of the guide-tube cannula and proximal sliding ol 

27 the second movable collar will extend the second extendable 

28 cannula or wire through the proximal end of the first 

29 extendable cannula. 

1 2. The device of claim 1, further comprising a 

2 first slide stop for the first movable collar that will limit 

3 the proximal distance of sliding of the first movable collar 

4 on at least one guide rail. 
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1 3. The device of claim 2, wherein the first slide 

2 stop is rigidly affixed to the first movable collar. 

1 4. The device of claim 2, wherein the first slide 

2 stop is rigidly affixed to at least one guide rail. 

1 5. The device of claim 1, wherein the second 

2 extendable cannula or wire is straight. 

1 6. The device of claim 1, wherein the second 

2 extendable cannula or wire is curved. 

1 7. The device of claim 1, further comprising a 

2 plate-like third or more movable collar mount slidably located 

3 on each guide rail rigidly affixed to a third extendable 

4 cannula having a proximal end and a distal end, wherein the 

5 distal end of the third extendable cannula is affixed to the 

6 third movable collar in a channel through the third movable 

7 collar and the proximal portion of the third extendable 

8 cannula is in the second extendable cannula; 

9 wherein proximal sliding of the third movable 

10 collar will extend the third extendable cannula through the 

11 proximal end of the second extendable cannula. 

1 8. The device of claim 1, wherein the first and 

2 second extendable cannulae comprise nitinol. 

1 9. The device of claim 1, further including a 

2 coordinated insertion means for pushing the foreign object 

3 into the body as it is pulled in by withdrawal of each wire or 

4 cannula in reverse order, said coordinated insertion means 

5 comprising: 

6 two concentric implant-bearing tubes of differing 

7 diameter encasing the foreign object, the larger of which is 

8 fastened to a block mounted on a stereotactic carrier such 

9 that th larger tub can be low r d onto th small r tube, 
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10 thereby advancing the implant under the guidanc of the 

11 smaller tube into an insertion trajectory in the body; 

12 a pulley-wheel fastened to the top of the 

13 stereotactic carrier on a rotationally adjustable mount, such 

14 that its axis is perpendicular to the line of the two 

15 concentric tubes and its plane of revolution can be adjusted 

16 to include the top of a second pulley-wheel on a rotationally 

17 adjustable mount fastened to the top of the device bearing the 

18 guide-tube cannula, curved first extendable cannula, and 

19 second extendable cannula or wire; 

20 a non-stretchable pulley cable fastened at one end 

21 to the block that drives the larger implant-bearing tube and 

22 extending to the pulley-wheel fastened to the top of the 

23 stereotactic carrier, further extending across to the pulley- 

24 wheel at the top of the device bearing the guide-tube cannula, 

25 curved first extendable cannula, and second extendable cannula 

26 or wire, and extending parallel to the axis of the guide-tube 

27 cannula to a fixation point in a second movable collar-mount 

28 which bears the second extendable cannula or wire; and 

29 a solid brace of adjustable length extending 

30 between the mounts of the rotationally adjustable pulley- 

31 wheels to maintain their positions when the non-stretchable 

32 pulley cable is pulled. 

1 10. A method for percutaneous placement of foreign 

2 substances in a desired location in a body, comprising: 

3 positioning a device as in claim 1 at a 

4 predetermined location and orientation over the body; 

5 advancing the guide-tube cannula into the body, 

6 covering an original predetermined distance; 

7 advancing the first extendable cannula into the 

8 body, covering a first additional predetermined distance; 

9 optionally advancing a second or more additional 

10 coaxial cannulae, covering a second or more additional 

11 predetermined distances; 

12 introducing the f reign substanc through th final 

13 ext ndable cannula t the desired location in the body. 
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1 11. The method of claim 10, wh rein the foreign 

2 substance is a catheter. 

1 12. The method of claim 10, wherein the foreign 

2 substance is a pharmaceutical compound. 

1 13. The method of claim 10, further comprising 

2 identifying the desired position in the body by an imaging 

3 technique. 

1 14. The method of claim 13, wherein the imaging 

2 technique is magnetic resonance imaging, computerized 

3 tomography, or ultrasound. 

1 15. A method for stereotactically implanting a 

2 foreign object along a compound trajectory through a desired 

3 region in a host, comprising: 

4 identifying an arcuate or helical curvilinear path 

5 through the host, which curvilinear path may traverse the 

6 region to be treated; 

7 inserting a guide-tube cannula through a surface of 

8 the host to meet the curvilinear path; 

9 advancing a curved cannula or wire, preformed to 

10 match the curvilinear path, through the guide-tube cannula 

11 into the host along the curvilinear path to a predetermined 

12 point; 

13 optionally advancing one or more additional coaxial 

14 cannulae, the last of which may be a solid wire, through the 

15 curved cannula to the surface of the host; 

16 attaching a foreign object to the final cannula or 

17 wire; 

18 withdrawing each wire or cannula through the host 

19 in reverse order, whereby the foreign object traverses the 

20 compound trajectory including the insertion path of each wire 

21 or cannula in the host. 
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1 16. The m thod of claim 15, wherein the foreign 

2 body is indirectly attached to th final cannula or wire. 

1 17. The method of claim 15, wherein the curved 

2 cannula or wire comprises nitinol. 

1 18. The method of claim 15, wherein the cannulae 

2 and wire are round, elliptical or rectangular in cross- 

3 section. 

1 19. The method of claim 15, wherein the compound 

2 trajectory is determined by a computer. 

1 20. The method of claim 15, wherein the foreign 

2 object includes a cryogenic probe. 

1 21. The method of claim 15, wherein the foreign 

2 object includes an electrode or electric device. 

1 22. The method of claim 15, wherein the foreign 

2 object is a catheter. 

1 23. The method of claim 22, wherein the catheter 

2 is semipermeable. 

1 24. The method of claim 23, wherein the catheter 

2 comprises polyacrylonitrile polyvinylchloride copolymer. 

1 25. The method of claim 23, wherein fluid 

2 withdrawn from the catheter is analyzed. 

1 26. The method of claim 23, wherein cells are 

2 present in the catheter. 
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1 27. Th method of claim 26, wherein the cells are 

2 c lis from the host that have be n transf ect d to express a 

3 therapeutic compound. 

1 28. The method of claim 23, wherein the catheter 

2 is f luidly connected to a source of therapeutic compound or 

3 biological agent. 

1 29. The method of claim 28, wherein the catheter 

2 is adjustably filled with the therapeutic compound or 

3 biological agent. 

1 30. The method of claim 29, wherein the disease is 

2 a tumor. 

1 31. The method of claim 29, wherein the disease is 

2 epilepsy, schizophrenia, Alzheimer's disease, depression, 

3 spinal cord trauma or Parkinson's disease. 

1 32. The method of claim 30, wherein the 

2 therapeutic compound is a dopamine agonist or a neurotrophic 

3 factor. 

1 33. The method of claim 22, wherein the catheter 

2 is a compound catheter, comprising: 

3 a central portion attached to connecting portions 

4 having tips inserted smoothly inside the central portion to a 

5 preselected point; 

6 a connector bead lodged in the inserted tip of each 

7 connecting portion and having a diameter large enough to press 

8 the outer wall of the connecting portion against the inner 

9 wall of the central portion, resulting in a tight, friction- 
10 fit seal; 

1X a terminal bead, friction-fit into one end of the 

12 compound cathet r, in which each end f said bead is 

13 internally shaped t a serrated c nical form, allowing a firm 
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14 friction-fit to a conplementarily shap d conical tip of a 

15 wire; and 

16 a wire passing through the length of the compound 

17 catheter and lodged in said terminal bead. 

1 34. A kit for placement of a catheter in a 

2 compound trajectory in a host, comprising: 

3 a device as in claim l; and 

4 a catheter. 

1 35. The kit of claim 34, wherein the catheter is 

2 semipermeable. 

1 36. The kit of claim 35, wherein the catheter 

2 comprises polyacrylonitrile polyvinylchloride copolymer. 

1 37. The kit of claim 34, further comprising at 

2 least one implantable pump. 

1 38. A method for determining a compound trajectory 

2 through a localised or extended region of a body for selective 

3 placement of a therapeutic modality, comprising: 

4 identifying the body region in a diagnostic image 

5 of tissue surrounding the region; 

6 assigning a coordinate system to the body; 

7 identifying three or more non-col linear points in 

8 the region; 

9 identifying the circle that intersects the three 

10 points, or a helix that passes near the four or more points; 

11 determining the radius, center and plane of the 

12 circle, or the formula of the helix; 

13 identifying points of intersection of the circle or 

14 helix with predetermined parallel sections through the body; 

15 identifying two points of tangency at which entry 

16 and exit segments of the compound trajectory contact an arc of 

17 the circle or helix; 
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18 determining the direction cosines of the entry and 

19 exit segments; 

20 calculating the arcuate distance between each pair 

21 of circle or helix intersection points on consecutive 

22 sections , from one point of tangency to the other; 

23 summing the arcuate distances; 

24 identifying points of intersection of the entry and 

25 exit segments with the predetermined parallel sections through 

26 the body; 

27 calculating the linear distances between each pair 

28 of entry or exit segment intersection points on consecutive 

29 sections; 

30 summing the linear distances; 

31 determining the compound trajectory of the path as 

32 the points of intersection of the arcuate and linear segments 

33 with the predetermined parallel sections; 

34 redefining the trajectory in the coordinate system 

35 of the diagnostic image; and 

36 displaying the trajectory in the diagnostic image. 

1 39. The method of claim 38, wherein the diagnostic 

2 image is a magnetic resonance image, computerized tomographic 

3 image or ultrasound image. 

1 40. The method of claim 38, wherein the body 

2 region is in the brain. 

1 41. The method of claim 40, wherein the body 

2 region is a tumor. 

1 42. The method of claim 40, wherein the body 

2 region is the basal ganglia. 
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